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ABSTRACT 

On the basis of the recently developed universal decline law of classical novae, we propose prediction for- 
mulae for supersoft X-ray on and off times, i.e., fx-on = (10 ± 1 .8)?3 days and tx-os = (5.3 ± 1 .4)(?3) 1 5 days for 
8 < f3 < 80 days. Here ?3 is the newly proposed "intrinsic" decay time during which the brightness drops by 
3 magnitudes from optical maximum along our universal decline law fitted with observation. We have deter- 
mined the absolute magnitude of our free-free emission model light curves and derived maximum magnitude 
vs. rate of decline (MMRD) relations. Our theoretical MMRD relations are governed by two parameters, 
one is the white dwarf (WD) mass and the other is the initial envelope mass at a nova outburst; this second 
parameter explains the scatter of MMRD points of individual novae. Our theoretical MMRD relations are 
also in good agreement with the well-known empirical formulae. We also show another empirical relation of 
My(15) ru — 5.7 ± 0.3 based on the absolute magnitude of our model light curves, i.e., the absolute magnitude 
at 15 days after optical maximum is almost common among various novae. We analyzed ten nova light curves, 
in which a supersoft X-ray phase was detected, and estimated their WD masses. The models best reproducing 
simultaneously the optical and supersoft X-ray observations are ONeMg WDs with 1.28 ±0.04 M (V598 
Pup), 1.23±O.O5M (V382 Vel), 1.15±O.O6M (V4743 Sgr), 1.13±0.06M© (V1281 Sco), 1.2±O.O5M 
(V597 Pup), 1.06 ±0.07 M (V1494 Aql), 1.04 ±0.07 M (V2467 Cyg), 1.07 ±0.07 M (V5116 Sgr), 
1 .05 ± 0.05 M (V574 Pup), and a CO WD with 0.93 ± 0.08 M (V458 Vul). The newly proposed relation- 
ships are consistent with the emergence or decay epoch of the supersoft X-ray phase of these ten novae. Finally, 
we discuss the mechanism of shock-origin hard X-ray component in relation to the emergence of companion 
star from the WD envelope. 

Subject headings: novae, cataclysmic variables — stars: mass loss — white dwarfs — X-rays: stars 



1. INTRODUCTION 

Classical novae are a thermonuclear runaway event on a 
white dwarf (WD) in a binary system, in which the WD ac- 
cretes hydrogen-rich matter from the companion star. When 
the accreted matter reaches a critical value, hydrogen at the 
bottom of the WD envelope ignites to trigger a shell flash 
(denoted by time ?ob)- Just after a nova outburst, the enve- 
lope on the WD rapidly expands to a giant size and optically 
thick winds blow. The WD envelope settles in a steady-state. 
The optical magnitude attains the peak at the maximum ex- 
pansion of the photosphere (time fo). Then the photosphere 
gradually shrinks whereas the total luminosity is almost con- 
stant during the outburst. Thus, the photospheric temperature 
r p h increases with time. The main emitting wavelength re- 
gion moves from optical to ultraviolet and finally to supersoft 
X-ray, corresponding to from 7J,h ~ 10 4 K through ~ 10 6 K. 
The optically thick winds blow continuously from the very 
beginning of the outburst until the photospheric temperature 
increases to log7ph (K)^ 5.4. Just when the optically thick 
wind stops (time f w ; n d), the temp erature quickly rises and the 
supersoft X-ray phase starts (e.g. JKato & Hachisu 111 9941) . 

Thus, classical novae become a transient supersoft X-ray 
source in a later phase of the outburst, but their X-ray detec- 
tion is rather rare ma inly because of sparse observing time of 
X-ray satellites (e.g jKrautter et al. Ifl996t lOrio et al. 112001 at 
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iNess et al. Il2007al) . If there is a formula to predict a supersoft 
X-ray phase of novae, we can plan more efficient observation. 
The aim of this paper is to propose such a formula. 

Novae show a rich variety of light curves; the timescales of 
decline, shapes, and total durations of outbursts are very dif- 
ferent from nova to nova. Therefore, it is difficult to predict 
a later supersoft X-ray phase only from the early optical light 
curve. Recently we have proposed the universal decline law of 
classical novae dHachisu & Kato II2006L 120071 lHachisu et al. I 
2008), i.e., there is a basic light curve which can basically fol- 
low (or overlap with) individual nova light curves by squeez- 
ing or stretching it in the direction of time. We call this "a 
universal decline law." The universal shape of light curves 
is almost independent of the WD mass and of the chemical 
composition of its hydrogen-rich envelope. This implies that 
each nova light curve can be uniquely specified by one param- 
eter (timescaling factor - squeezing/stretching factor against 
a universal shape of nova light curves). 

We have reproduced a number of light curves of novae 
from e arly phase to supersoft X-ray phase (turn-on/turnof f 
times) (Ha chisu&Katol 120061 12001; lHachisu et al71 120061 
120071 120081) . Based on these experiences we are now able to 
construct empirical formulae that tell us the emergence/decay 
times of a supersoft X-ray phase of individual novae from ob- 
servational information of the early phase. We hope that these 
formulae are useful for future X-ray observations of classical 
novae. 

The next section (Section [2]) describes our numerical 
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method for obtaining X-ray, ultraviolet (UV), optical, and in- 
frared (IR) light curves. In Section [3] we show homologous 
shapes of our model light curves. Using characteristic proper- 
ties of free-free emission model light curves, we derive the ab- 
solute magnitude of each optical light curve, which provides a 
unique standard candle of novae. We further propose another 
method for obtaining the nova absolute brightness by fitting 
a nova light curve with the other the brightness of which is 
known. After various nova timescales and relations between 
these timescales are presented, we propose our prediction for- 
mulae of a supersoft X-ray phase in Section [4] Then, we 
present our theoretical understanding on the so-called max- 
imum magnitude vs. rate of decline (MMRD) relation of 
classical novae in Section [5] The light-curve analyses of ten 
classical novae are presented in Section [6] The X-ray turn- 
on/turnoff times of these ten novae, i.e., V598 Pup, V382 Vel, 
V4743 Sgr, V1281 Sco, V597 Pup, V1494 Aql, V2467 Cyg, 
V51 16 Sgr, V574 Pup, and V458 Vul, are consistent with the 
above relations (formulae). Conclusions follow in Sections[7] 

2. MODELING OF NOVA LIGHT CURVES 

The decay phase of nova e can be followed by a s equence of 
steady state solutions (e.g jKato & Hachisu 119941). Using the 
same m ethod and numerical techniques as in lKato & Hachisu I 
GUI, we have followed evolutions of novae by connecting 
steady state solutions along the decreasing envelope mass se- 
quence. The mass of the hydrogen-rich envelope is decreasing 
due to wind mass-loss and nuclear burning. We solve a set of 
equations, consisting of the continuity, equation of motion, 
radiative diffusion, and conservation of energy, from the bot- 
tom of the hydrogen-rich envelope through the photosphere 
assuming spherical symmetry. Winds are accelerated deep in- 
side the photosphere so that they are called "optically thick 
winds." 

We have presented a unified model for IR, optical, UV, 
and super soft X-ray light curves of the decay phase of classi- 
cal novae (lHachisu & Kato II2006L 1200 71: lHachisu et al. II2007L 
I2008t iKato & Hachisu 1 120071) . The optical and IR luminos- 
ity are repro duced by free-free emiss i on from optically thi n 
plasma (e.g., iGallagher & Nev Ifl976t iKrautter et al. 1 1 19841) . 
This free-free emiss ion flux at the frequency v is estimated 
from Equation (9) of lHachisu & Kato I (120061) . i.e., 
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where M w i n d is the wind mass-loss rate, v p h the wind veloc- 
ity at the photosphere, and R p h the photospheric radius, all of 
which are taken from our wind solutions. It should be noted 
that the shapes of free-free emission light curves are indepen- 
dent of the frequency (wavelength). After the optically thick 
wind stops, the total mass of the ejecta remains constant in 
time. The flux from such hom ologously expanding eje cta is 
estimated from Equation (19) of Hachis u & Kato I (120061) . i.e., 



F v oc t 



(2) 



where t is the time after the outburst. 

Assuming a blackbody spectrum with the photospheric 
temperature of 7ph, we have calculated the UV 1455 A light 
curv e for a narrow energy band of 1445-1465 A (see, 
e.g., ICassateUa et"aTll2002]; lHachisu & Kato 1120061 for more 
details) and have also estimated the supersoft X-ray flux 
dKato & Hachisu I [19941 lHachisu & Kato 1 120061) throughout 



TABLE 1 

Chemical Composition of White Dwarf 
Envelope 



nova case 


X 


Y 


XcNO 




Z a 


CO nova 1 


0.35 


0.13 


0.50 


0.0 


0.02 


CO nova 2 


0.35 


0.33 


0.30 


0.0 


0.02 


CO nova 3 


0.45 


0.18 


0.35 


0.0 


0.02 


CO nova 4 


0.55 


0.23 


0.20 


0.0 


0.02 


CO nova 5 


0.65 


0.27 


0.06 


0.0 


0.02 


Ne nova 1 


0.35 


0.33 


0.20 


0.10 


0.02 


Ne nova 2 


0.55 


0.30 


0.10 


0.03 


0.02 


Ne nova 3 


0.65 


0.27 


0.03 


0.03 


0.02 


Solar 


0.70 


0.28 


0.0 


0.0 


0.02 


a carbon, nitrogen, oxygen, and 


neon are 


also included in 



Z = 0.02 with the same ratio as the solar abundance 

the present paper for t he energy range of 0.2-0.6 keV (e.g., 
lHachisu & Kato ll2009l) . Even if we adopt a different energy 
range of 0.2-2.0 keV, shapes of supersoft X-ray light curves 
hardly change because the blackbody photosphere emits pho- 
tons with energy of < 0.6 ke V for Mwd ^1-2 M Q . 

The light curves of our optically thick wind model are pa- 
rameterized by the WD mass (Mwd), chemical composition of 
the WD envelope {X, Y, Xcno, ^Ne, and Z), and envelope mass 
(M env o) at the time of the outburst (?ob)- Here X is the hydro- 
gen, Y the helium, Xcno the carbon-nitrogen-oxygen, X^ e the 
neon, and Z the metallicity content of the envelope. We adopt 
Z = 0.02, which also includes carbon, nitrogen, oxygen, and 
neon with solar composition ratios. Table Q] shows nine sets 
of chemical compositions assumed in our models. The details 
of our mode l light curves and t heir p arameter dependence are 
described in Hachis u & Kato I (120061) . 

3. UNIVERSAL DECLINE LAW 

The decay timescale of a light curve depends mainly on the 
WD mass and weakly on the chemical composition, whereas 
the maximum brightness of a light curve depends strongly on 
the initial envelope mass, M env q. In this section, we present 
various features, especially the absolute magnitudes, of nova 
light curves based on our universal decline law. 

3.1. Template of nova light curves 

Figure Q] shows our free-free emission model light curves 
calculated from Equation (Q]), i.e., 

-| {*?wd} 

+ G {Mwd} , (3) 



nig = -2.5 log 



M wind 



1 ph 



Rph 



for various WD masses, from Mwd = 0.55 to 1 .2 M Q by step 
of 0.05 M for a chemical composition of "CO nova 2" (these 
numerical data are tabulated in Table|2|. The subscript (f) ex- 
plicitly expresses that this is a function of time while the su- 
perscript {Mwd} indicates a model parameter of WD mass. 
An optically thick wind phase continues for a long time, from 
a very early phase to a late phase of a nova outburst. In this 
figure we plot free-free emission model light curves calcu- 
lated by Equation (0 from a very early phase until the end of 
an optically thick wind phase. The 15th mag point of each 
light curve corresponds to the end of the wind phase. In other 
words, we define a constant of G^ MwD ^ in Equation (O as that 
the last (lowest) point of each light curve (the end of an opti- 
cally thick wind phase) is 15th mag. As a result, the value of 
^{Mwd} j s different for different Mwp- In order to make Table 
|2]more compact, we have defined G^ MwD ^ in this way. 
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TABLE 2 
Light Curves of CO Novae' 





0.55M Q 


O.6M 


0.65M Q 


O.7M 


O.75M 


O.8M 


O.85M 


O.9M 


O.95M 


l.OAf 


l.O5M 


1.1M 


1.15M 


\.1Mq 


(mag) 


(day) 


(day) 


(day) 


(day) 


(day) 


(day) 


(day) 


(day) 


(day) 


(day) 


(day) 


(day) 


(day) 


(day) 


1.50 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


1.75 


1.319 


1.421 


1.286 


1.319 


1.231 


1.187 


1.123 


1.120 


1.032 


0.970 


0.922 


0.870 


0.772 


0.717 


2.00 


2.921 


2.835 


2.615 


2.705 


2.442 


2.367 


2.243 


2.210 


2.027 


1.910 


1.811 


1.660 


1.528 


1.428 


2.25 


4.634 


4.393 


4.162 


4.102 


3.737 


3.548 


3.343 


3.220 


2.997 


2.850 


2.660 


2.440 


2.247 


2.140 


2.50 


6.346 


6.132 


5.758 


5.482 


5.069 


4.750 


4.453 


4.220 


3.917 


3.780 


3.485 


3.200 


2.940 


2.820 


2.75 


8 187 


7 872 


7 358 


6 842 


6 398 


5 960 


5 543 


5 220 


4 837 


4 700 


4 280 


3 940 


3 622 


3 470 


3.00 


10.46 


9.822 


9.138 


8.272 


7.738 


7.187 


6.623 


6.220 


5.747 


5.620 


5.057 


4.650 


4.276 


4.100 


3.25 


12.72 


12.05 


1 1.06 


9.992 


9.088 


8.409 


7.713 


7.170 


6.627 


6.480 


5.797 


5.340 


4.926 


4.710 


3.50 


16.25 


14.32 


12.99 


1 1.71 


10.55 


9.659 


8.793 


8.130 


7.517 


7.260 


6.527 


6.010 


5.548 


5.310 


3.75 


21 22 


17.04 


15 39 


13 50 


12 15 


10 98 


9 913 


9 100 


8 407 


8 060 


7 257 


6 670 


6 169 


5 890 


4.00 


27 67 


19.96 


17 88 


15.66 


13 80 


12 32 


1 1 05 


10 07 


9 257 


8 850 


7 967 


7 320 


6 774 


6 470 


4.25 


34 30 


23 75 


20 45 


17 95 


15 47 


13.69 


12 21 


1 1.06 


10 13 


9 650 


8 677 


7 970 


7 369 


7 040 


4.50 


41.31 


28.68 


23.06 


20.29 


17.18 


15.20 


13.42 


12.08 


1 1.02 


10.44 


9.387 


8.620 


7.969 


7.610 


4.75 


48.32 


34.46 


27.56 


22.96 


19.69 


16.79 


14.77 


13.13 


1 1.97 


11.24 


10.11 


9.280 


8.579 


8.160 


5 00 


55 67 


42 32 


33 29 


26.67 


22 73 


18 44 


16 29 


14.27 


12 95 


12 09 


10 87 


9 950 


9 209 


8 740 


5.25 


65.40 


50.31 


39.71 


32.68 


25.92 


21.12 


17.85 


15.78 


14.08 


12.99 


1 1.67 


10.66 


9.899 


9.340 


5.50 


75.12 


58.89 


46.20 


38.40 


29.35 


24.37 


20.06 


17.41 


15.58 


13.98 


12.56 


11.44 


10.73 


10. 1 1 


5.75 


87.91 


68.64 


52.69 


43.44 


33.25 


27.64 


23.20 


19.35 


17.32 


15.21 


13.65 


12.51 


11.60 


10.99 


6.00 


101 8 


78 86 


60 29 


48 43 


37 47 


30 94 


25 89 


21.74 


19 04 


16.97 


15 15 


13 97 


12 77 


1 1 93 


6.25 


1 17.5 


90.97 


68.44 


54.61 


41.95 


33.84 


28.14 


23.89 


20.85 


18.59 


16.54 


14.98 


13.83 


12.90 


6.50 


134.6 


103.7 


76.86 


61.70 


46.86 


37.24 


30.55 


25.75 


22.50 


19.92 


17.56 


15.92 


14.53 


13.59 


6.75 


152.9 


1 17.3 


87.51 


69.39 


52.1 1 


41.14 


33.37 


27.75 


24.15 


21.28 


18.61 


16.71 


15.21 


14.14 


7.00 


171.3 


131.5 


98.25 


77.79 


57.84 


45.35 


36.66 


30.1 1 


25.92 


22.68 


19.72 


17.57 


15.94 


14.7 1 


7.25 


190.6 


146.1 


109.6 


86.81 


64.22 


49.91 


40.23 


32.84 


27.97 


24.18 


20.91 


18.58 


16.74 


15.34 


7.50 


211.6 


161.1 


121.3 


96.18 


71.09 


54.93 


44.08 


35.79 


30.29 


25.87 


22.23 


19.67 


17.58 


16.03 


7.75 


234.6 


177.5 


133.6 


105.6 


78.50 


60.32 


48.19 


38.77 


32.77 


27.76 


23.70 


20.82 


18.49 


16.75 


8.00 


260 1 


196 


146.7 


115 7 


86 05 


66 19 


52 32 


41 90 


35 42 


29 80 


25 23 


22 03 


19 46 


17 49 


8.25 


288 9 


217 2 


160.5 


127 1 


94 01 


72 60 


56 72 


45.29 


38 05 


31 91 


26 86 


23 28 


20 44 


18 30 


8 50 


322 8 


241 2 


177 7 


139 7 


102 9 


79 60 


61 60 


48 95 


40 89 


34 14 


28 56 


24 53 


21 47 


19 12 


8.75 


363.7 


268.6 


197.2 


153.9 


113.1 


87.22 


66.89 


53.02 


44.02 


36.54 


30.38 


25.89 


22.58 


19.97 


9.00 


411.1 


301.4 


219.4 


171.6 


124.3 


95.52 


73.07 


57.43 


47.44 


39.14 


32.31 


27.46 


23.81 


20.94 


9.25 


457.6 


339 


246 2 


191 7 


138 9 


104 5 


80 58 


62 78 


51 34 


42 10 


34 43 


29 22 


25 19 


22 04 


9.50 


508.2 


377.0 


276.1 


214.1 


155.7 


1 17.0 


89.01 


69.27 


56.25 


45.37 


37.10 


31.32 


26.89 


23.32 


9 75 


547 2 


412 


307 1 


239 8 


173 7 


131 4 


100 


76 87 


62 24 


49 59 


40 33 


33 85 


28 93 


24 79 


10 00 


589 5 


446.5 


335.6 


267 7 


193 2 


146 7 


112 5 


86.46 


69 49 


54 71 


44 21 


37 08 


31 56 


26.67 


10.25 


635.0 


479.3 


362.2 


292.3 


215.6 


163.5 


124.7 


97.02 


78.1 1 


60.81 


48.97 


40.76 


34.69 


29.01 


10.50 


673.6 


512.3 


386.2 


313.5 


234.6 


180.4 


138.3 


107.1 


86.93 


68.04 


54.70 


45.39 


38.36 


31.74 


10.75 


715.8 


541.6 


412.6 


335.6 


253.3 


195.2 


153.4 


1 18.5 


95.91 


75.64 


60.98 


50.50 


42.51 


34.85 


1 1 .00 


762.2 


571.6 


440.2 


355.9 


269.9 


209.3 


165.7 


130.3 


105.7 


83.14 


66.78 


55.83 


46.53 


38.24 


11.25 


813.6 


603.7 


464.6 


377.3 


284.6 


223.5 


177.7 


140.4 


114.3 


91.54 


73.13 


60.75 


50.34 


41.48 


11.50 


869.2 


637.0 


489.4 


396.0 


300.8 


237.0 


187.3 


149.9 


121.1 


99.24 


79.95 


65.55 


54.40 


44.25 


11.75 


920.8 


673.7 


516.6 


415.2 


317.0 


249.3 


197.7 


158.9 


128.5 


105.0 


84.82 


70.35 


58.73 


47.19 


12.00 


975.4 


713.8 


546.4 


436.2 


332.6 


261.8 


209.0 


167.0 


136.5 


111.3 


89.90 


74.95 


62.18 


50.27 


12.25 


1034. 


756.3 


578.9 


459.2 


349.7 


274.1 


219.3 


175.8 


145.2 


118.1 


95.36 


79.25 


65.64 


53.55 


12.50 


1095. 


801.4 


613.5 


484.3 


368.4 


287.4 


229.3 


184.0 


154.3 


125.4 


101.1 


83.80 


69.32 


56.94 


12.75 


1160. 


849.1 


650.2 


513.2 


388.7 


302.1 


240.2 


192.7 


161.7 


133.1 


107.3 


88.40 


73.16 


60.21 


13.00 


1229. 


899.8 


689.0 


544.1 


411.5 


317.9 


252.1 


202.1 


169.7 


139.7 


113.6 


93.20 


77.25 


63.71 


13.25 


1302. 


953.3 


730.1 


576.7 


436.1 


336.6 


265.4 


212.3 


178.3 


146.4 


120.7 


98.30 


81.99 


67.40 


13.50 


1379. 


1011. 


773.7 


611.3 


462.2 


356.7 


281.6 


224.3 


187.5 


153.8 


128.1 


103.8 


87.01 


71.28 


13.75 


1461. 


1071. 


819.8 


647.9 


489.8 


377.9 


298.9 


238.4 


199.1 


161.6 


136.0 


109.6 


92.33 


75.31 


14.00 


1547. 


1134. 


868.7 


686.7 


519.0 


400.5 


317.1 


253.2 


211.4 


173.1 


144.4 


116.8 


97.94 


80.17 


14.25 


1639. 


1202. 


920.5 


727.8 


550.0 


424.3 


336.5 


269.0 


224.4 


185.5 


153.2 


124.8 


103.9 


85.41 


14.50 


1736. 


1273. 


975.3 


771.3 


582.9 


449.6 


357.0 


285.7 


238.2 


198.6 


162.6 


133.3 


110.2 


90.96 


14.75 


1839. 


1349. 


1034. 


817.4 


617.7 


476.4 


378.6 


303.4 


252.9 


212.4 


172.6 


142.3 


116.9 


96.87 


15.00 


1948. 


1429. 


1095. 


866.2 


654.5 


504.8 


401.6 


322.1 


268.3 


227.1 
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a chemical composition of the envelope is assumed to be that of "CO nova 2" in TableQ 
k duration of supersoft X-ray phase in units of days 
c converted magnitudes at the bottom point in Figure [2l 

stretching factors against VI 668 Cyg observational data in Figure[2l 
c absolute magnitudes at the bottom point in Figure [2|by assuming (m— M)y = 14.3 (V1668 Cyg) 
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TABLE 3 
Light Curves of Ne Novae 2 
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a chemical composition of the envelope is assumed to be that of "Ne nova 2" in Table ITI 
k duration of supersoft X-ray phase in units of days 
c converted magnitudes at the bottom point in Figure [6| 

stretching factor against V1974 Cyg observational data in Figure [6| 
c absolute magnitudes at the bottom point in Figure[6]by assuming (m—M)v = 12.3 (V1974 Cyg) 
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days after outburst (+ const.) 

FIG . 1 . — Magnitudes of our free-free emission model light curves for various WD masses, i.e., 0.55 — 1.2 Mq by 0.05 Mq step, numerical data of which are 
tabulated in Table[2] We adopt a chemical composition of "CO nova 2" in Table [l] for these WD envelopes. The decay timescale depen ds mainly on th e WD 
mass. We add observational y magnitudes (open squares) of V1668 Cyg (Nova Cygni 1978) for comparison, which are taken from Gallagher et al. 1 1980). The 
y magnitudes are shifted up by 1.2 mag to fit them with the 0.95 Mq WD model (see the right axis). Our light curve model of 0.95 Mq WD (blue color line) 
almost perfectly fits with the observational y magnitudes except the last point. 
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t'=t/f s (days) 

FIG. 2. — Same as Figure[T] but in the plane of rescaled flux, —2.5 log/sF^ +const. =my —2.5 log / s =mL and scaled time, t/f s =t'. We added UV 1455A model 
light curves which are stretched along time by a factor of l//s. This stretchin g factor is determined in a way that each UV light curve is overlapped to that of 
V1668 Cyg observation (red big open circles), data of which are taken from Cassatella et al. (2002). The UV model fluxes are normalized to the peak value 
of the UV observation. The values of log/ s are tabulated in Table [2] The sudden drop in the UV flux of VI 668 Cyg on day ~ 60 is caused by formation of 
an optically thin dust shell (denoted by an arrow with "dust"). All the free-free emission model light curves are almost overlapped to each other. We call this 
homologous nature of free-free light curves "the universal decline law." Here t% and times are not uniquely determined along our universal decline law, because 
these times depend on the initial envelope mass as indicated by two sets of ((o, ti, (3) in the figure. Three different initial envelope masses of M em q = 2.3 X 10~ 5 , 
1.6 X 10~ 5 , and 1.1 X 1O~ 5 M0 correspond to three peak brightnesses, point s A, B, and C, respectively, o n the 0.95 Mq WD model. We added the V magnitude 
observation (red filled circles) for comparison, data of which are taken from Mallama & Skillman 1 1979). The y magnitudes are shifted up by 0.3 mag to match 
them with the V magnitudes in the early phase. 
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These light curves demonstrate that the decay timescales 
of novae are sensitive to the WD mass. Therefore we can 
determine the WD mass of a nova from light curve fitting. 
The 0.95 Mq (blue thick solid line) or 1 .0 Mq WD may be a 
typical case of fast novae and the light curve of 0.65 Mq WD 
(red thick solid line) may correspond to a slow nova. 

These free-free emission model light curves are applica- 
ble to optical and near i nfrared (IR) wavelength regions (e.g., 
lHachis u & Katol l2006T) . We added observational y magni- 
tudes of V1668 Cyg (Nova Cygni 1978) to Figure Q] which 
nicely reproduce the light curve of 0.95 Mq WD. In this case, 
we shift the model light curve of 0.95 Mq WD down by 1.2 
mag in order to fit it with the observation (see the right axis 
of Figure [TJ. The intermediate width y band is an emission 
line-free optical band, so it is ideal to fit our model light curve 
with y magnitudes. 

In this way, we can shift the template light curves to fit them 
with observation for individual novae. In other words, one can 
shift the model light curve "up and down" against the obser- 
vational data of an individual nova. In the case of V1668 Cyg, 
we obtain m y = ma+ 1 .2 from the 0.95 Mq WD model in Table 

m 

3.2. Brightness of free-free light curves 

The free-free emission model light curves in Figure [T]have 
strong similarity in the ir shapes. In our previous paper 
(lHachisu & Kato 1 12006). we showed that these model light 
curves are homologous and their shapes are almost overlapped 
with each other when they are properly squeezed or stretched 
along time. In this subsection, we reformulate this homol- 
ogous properties in a more sophisticated manner and obtain 
their brightnesses. 

First, we squeeze or stretch our UV 1455 A model light 
curves in the direction of time and also its peak flux is normal- 
ized to be fitted with the observational peak of V1668 Cyg. 
Then we obtain all the UV light curves almost overlapped to 
each other as shown in Figure [2] We determine a stretch- 
ing factor, / s , of each UV model light curve by increasing or 
decreasing it until the model light curve shape matches the 
observation. Table [2] shows logarithmic values of log/ s , that 
is, positive (negative) for stretching (squeezing) against the 
V1668 Cyg UV 1455 A observation. For example, we see 
that the 0.95 Mq WD model light curve is almost overlapped 
to the VI 668 Cyg observation without stretching (/ s w 1.0) 
and the 0.8 Mq WD model light curve evolves about twice 
slower (/ s w 2.0) than the 0.95 Mq WD model. In this way, 
all the UV 1455 A model light curves are overlapped to each 
other, if we squeeze/stretch the time as t' = t / f s . Note that 
the stretching factor is slightl y different from th at in our pre- 
vious paper (Table 9 of Hachisu & Kato 1 12006), which was 
determined in the way that each free-free model light curve 
is overlapped with each other instead of UV 1455 A model 
light curves. 

Figure [2] also shows optical data (y and V magnitudes) of 
VI 668 Cyg and fitted free-free emission model light curve of 
0.95 Mq WD (/ s w 1.0). The other free-free emission model 
light curves are squeezed/stretched with the same / s as that of 
UV light curves. When we squeeze the model light curves by 
a factor of / s in the time direction (f ' = t // s and v' = f s v), the 
flux is also changed as F' v , = f % F v because 



Therefore, we shift the model light curves in the vertical di- 
rection as 



d_ = d 

dt' **** 



■ m v - 2.5 log/ s , 



(5) 



and obtained Figure [2] Here my is the apparent V magni- 
tude of an Mwd model and m' v is its time-stretched magnitude 
against a standard model of 0.95 M© WD (« V1668 Cyg). 
We see that all the light curves are almost overlapped to each 
other. For more details of derivation of Equation (O, see Ap- 
pendix A. 

The origins of magnitude in Figures Q] and [2] are different 
because the magnitude in the right axis of Figure Q] is cali- 
brated by the y magnitudes whereas by the V magnitudes in 
Figure|2] We define the rescaled apparent V magnitude (time- 
stretched V magnitude) against the V magnitudes of VI 668 
Cyg instead of the y magnitudes. The apparent V magnitudes 
are written as 



m' v = -2.5 log 



-i {M W d} 



+ K V 



(6) 



J (t/f.) 



where Ky is a constant and is determined to match m' v of the 
0.95 M Q WD model with the V magnitude of V1668 Cyg as 
shown in Figure|2] Therefore, Ky is also related to G^ 95 

as 



This simply means 



= my = mff+0.9, 



(8) 



dt 



(4) 



for 0.95 M WD model with / s = 1.0. Note that we here ob- 
tain 0.9 instead of 1.2 because observational m v is brighter 
than to v by about 0.3 mag, that is, my and m y -Q3 are over- 
lapped to each other in Figure [2] We use the same Ky for all 
other WD mass models (fromM W D = 0.55 to 1.2 M by 0.05 
Mq step) in Figure [2] All the free-free emission model light 
curves are almost overlapped to each other. 

It should be noticed here that the flux is condensed and in- 
creased by a factor of / s when the time is stretched by a factor 
of 1 /f s (or squeezed by a factor of / s ) only under the condition 
that the spectrum is not changed after time-stretching. This 
condition is satisfied in free-free emission (optical and near 
IR) because the flux (F u ~ constant) is independent of the fre- 
quency v, but not satisfied in blackbody emission because the 
blackbody emissivity depends on the frequency v and, after 
time-stretching, a blackbody spectrum changes with v' = f s v. 
See Appendix A for more detailed formulation of the scaling 
of free-free flux. 

3.3. Maximum brightness vs. initial envelope mass 

Figure|2]shows three points of A, B, and C on the same light 
curve of 0.95 Mq WD model. Point B indicates the initial 
envelope mass of M env o = 1.6 x 10~ 5 M Q , corresponding to fo 
at the maximum brightness of V1668 Cyg. Point A indicates 
a much brighter maximum and corresponds to a much larger 
envelope mass (M env .o = 2.3 x 10~ 5 Mq) than that of point B. 
On the other hand, point C corresponds to a much smaller 
envelope mass (M env ,o = 1.1 X 10~ 5 Mq). 

The initial envelope mass depends on the mass accretion 
rate to the WD before ignition ( see, e.g., Figure 9 of Nomoto 
Il982t iPrialnik & Kovetz1ll995h . Point A is corresponding to 
a case of much lower mass accretion rates to the WD. On 
the other hand, point C is corresponding to a case of much 
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FIG. 3. — Same as Figure[2] but in the logarithmic time vs. magnitude plane. The origin of time is about 4 days before the optical maximum of V1668 Cyg. 
The right edges of free-free emission model light curves correspond to the epoch when the optically thick winds stop (denote by an arrow labeled f w i„d). We also 
added a decline rate of free-free flux, i.e., <x f - ' 75 (blue dash-dotted line) in the early phase before fbreak- 



Absolute Magnitude 




12 3 

log t (days) 

FIG. 4. — Same as Figure[3] but for absolute magnitudes and real timescales. We have calibrated the free-free model light curves by a distance modulus of 
(m—M)v = 14.3 for V1668 Cyg. The position at point B in Figure[3]is indicated by a blue filled circle. We also show the magnitude, My(15), at 15 days after the 
optical maximum by a blue cross. Their average value of My(15) = -5.95 ± 0.25 is also indicated in the figure for 0.7-1.05 Mq WDs. Model UV and supersoft 
X-ray light curves are omitted to simplify the figure except three WD masses (blue solid line for 1.0, red solid line for 0.95, and green solid line for 0.9 Mq). 
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higher mass accretion rates to the WD than that of point B, 
for the case of V1668 Cyg. This simply suggests that novae 
are brighter for lower mass accretion rates. 

For the other WD masses, the trend between the maximum 
brightness and the initial envelope mass is the same as that of 
0.95 Mq. For more (less) massive WDs, the envelope mass at 
the same point B is smaller (larger) than that of 0.95 Mq. We 
will see again these relations between the maximum bright- 
ness and the initial envelope mass in more details in Section 

s 

We regard that VI 668 Cyg is a typical classical nova and 
point B in Figure |2] corresponds to typical brightnesses for 
other WD mass models. Therefore, we replot the same nova 
model light curves in a logarithmic time in Figure [3] In the 
overlapping procedure of UV 1455A model light curves, we 
adopt the origin of time which is different from the origin of 
time in raff of Table [2] This procedure is necessary for the 
horizontal axis of log?' because we must start all the free-free 
emission light curves from the same stage like point B (or 
C) in Figure [2] Here we define the origin of time as about 4 
days before point B in Figure [2] that is, the outburst day of 
V1668 Cyg. In this case, the initial brightnesses of novae are 
almost the same (brightness at point B) and overlapped with 
each other in the \ogt'-m' v diagram. 

3.4. Absolute magnitude of free-free light curves 

We have given an apparent magnitude to each free-free 
emission model light curve against VI 668 Cyg data. Next 
step is to determine the absolute magnitudes of these model 
light curves. This can be done by using the distance modulus 
of V1668 Cyg. Figure [2] shows observed y and V magnitudes 
of VI 668 Cyg, which can be converted to the absolute mag- 
nitude from the distance modulus to V1668 Cyg, i.e., 



(m-M) v = [51og(rf/10)+A v 



I V 1668 Cyg 



= 14.3. 



(9) 



Here, we adopt the distance of d = 4. 1 kpc dKato & Hachisu I 
120071) and the absorption of Ay = 3.1 x E(B-V) = 3.1 x 
0.4 = 1.24 dSfickland et al. Ifl98ll) for V1668 Cyg. Since our 
0.95 Mq WD model shows a good agreement with the y and 
V magnitudes of VI 668 Cyg, we determine the absolute V 
magnitude of our 0.95 M© WD model light curve as follows: 



My = Illy — 14.3 = (ffIff+0.9)- 14.3, 



(10) 



where m v = m v = ra ff + 0.9 for 0.95 Mq WD model (/ s = 1) 
from Equation (0. Then, the absolute V magnitude at the last 
point of free-free template light curve (the end point of Table 
13 corresponds to 



M w = m w - 14.3 = (ra ff + 0.9) -14.3= 1.6, 



(11) 



where the subscript of "w" means "wind" which corresponds 
to rag = 15 at the end point (Tabled- Using M w , we restore 
the entire absolute V magnitudes of the 0.95 M© WD model 
by 

My = (ra ff - 15) +M W , (12) 

and it is shown by a red solid line in Figure H] 

For the other WD mass models, from Equation (0 we sim- 
ply derive the following relations, i.e., 



My = ra v + 2.5 log/ s - 14.3, 
M w = ra; + 2.51og/ s -14.3, 



(13) 
(14) 
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FIG. 5. — The absolute magnitude at the end point of each free-free emis- 
sion model light curve, M w , is plotted against the WD mass for two differ- 
ent chemical compositions, i.e., "CO nova 2" (star marks) and "Ne nova 2" 
(squares). See text for details. 

based on the VI 668 Cyg data. We plot M w against various 
WD masses in Figure|5] Using M w , we have 



My = (ra ff - 15) + M W , 



(15) 



where log/ s , ra^, and M w are tabulated in Table [2] for all 
other WD mass models from Mwd = 0.55 to 1.2 Mq by 
0.05 Mq step. Here prime symbol means the values for 
squeezed/stretched light curves. We have restored the entire 
absolute V magnitudes for 0.55-1 .2 Mq WDs in Figure|4] 

It is interesting to examine whether or not the absolute 
magnitudes of our free-free emission model light curves are 
consistent with various empirical relations of classical no- 
vae that have ever been proposed. Very popular relations 
between the maximum magnitude and the rate of decline 
(MMRD) will be discussed in Section [5] Here we check 
the empirical formula that the absolute magnitude at 15 days 
after optical maximum , My(15), is almost common among 
various novae (e.g., iBuscombe & de Vaucouleurs I 119551 



I Capaccioli et al. Nl989UCohenl 119851: iDownes & Duerbeckl 



2000; van den Bergh & Youn ger 1 119871) 
first proposed by Buscombe & de Vau couleurs 
M v (15) = -5.2 ± 0.1, the n followed bvlCohen 



This relation was 
1955]) with 
1985) with 



M v (15) = -5.60 ± 0.43, bv lvan den Bergh & Younger I (1 19871) 
with M V (15) = -5.23 ±0.39. bv lCapaccioli et al. I d!989b with 
M V (15) = -5.69 ± 0.42, and bv lPownes & Duerbeckl ( 120001) 
with M V (15) = -6.05 ± 0.44. We have obtained M y (15) = 
-5.7 ±0.3 for 0.55-1. 2 M Q WDs (14 light curves) with equal 
weight for all WD mass models, being roughly consistent 
with the above empirical estimates as shown by blue crosses 
in Figure 2] For the statistical point of view, carbon-oxygen 
(CO) novae have typical masses of 0.7-1.05 Mq. If we take 
this rage of WD masses (equal weight), the value becomes 
My(15) = -5.95 ±0.25 a s indicated in Figur e |H being close 
to the value proposed bv lPownes & Duerbeck (200(|. 



3.5. Absolute magnitudes for neon novae 

In this subsection, we study light curves of neon (Ne) novae 
and show that the essential results on the free-free emission 
model light curves are the same as those for CO novae studied 
in the previous subsections. Table[3]lists free-free model light 
curves in the form of raff, definition of which is the same (see 
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log(t/f s ) (days) 

FIG. 6. — Same as Figure[3] but models for "Ne nova 2." Observational data of the neon nova V1974 Cyg is added for comparison. Thirteen different WD 
mass models (0.7- 1.3 Mq by 0.05 Mq step) are plotted in the logarithmic time-magnitude plane, data of which are tabulated in Table[3] We squeeze or stretch 
the model light curve along time to fit with the observational data of V1974 Cyg, i.e., ROSAT X-ray observation, optical, and IUE UV 1455A light curves. UV 
1455A model fluxes are normalized to have the same peak value as that of IUE data, peak of which is 1.5 4 X 10~ n erg s~ * cm -2 A -1 . Fluxes of supersoft X-ray 
are normalized to fit it with the peak of ROSATX-my data. The observational V magnitudes are taken from Chochol et al. ( 1993), visual data from the American 
Association of Variable Star Observers (AAVSO), IUE UV 1455 A data from lCassatella et al. I <2"00"2T) . and ROSAT X-ray data from lKrautter et al. I fT99H) . We 
also added two decline rates of free-free flux, i.e., Fx oc t (blue dash-dotted line) in the early phase before t^ rs .. d ^, and Fx oc t' 3 (blue solid straight line) of 
Equation (2} in the expanding nebular phase. Here we set point B on each model light curve near tg . 



t 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 — I 6 




log t (days) 

FIG. 7. — Same as Figure|6] but for absolute magnitudes and real timescales. We have calibrated the free-free model light curves by a distance modulus of 
(m-M)v = 12.3 for V1974 Cyg. The positions at point B in Figure[6]are indicated by a magenta filled circle. We also show the magnitude, Mv(15), at 15 days 
after the optical maximum by a magenta cross. Their average value of My (15) = —5.85 ± 0.09 is also indicated in the figure for 0.9-1.15 Mq WDs. Model UV 
and supersoft X-ray light curves are omitted to simplify the figure except three WD masses {blue solid line for 1.1, red solid line for 1 .05, and green solid line for 
1.0 Mq). 
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Equation (O) but for the chemical composition of "Ne nova 
2" in Table Q] 

Figure [6] shows free-free (optical), UV 1455 A, and super- 
soft X-ray model light curves as well as observational data of 
V1974 Cyg (Nova Cygni 1992). We plot these light curves 
in log t '-m' v plane and show nova evolutions until more later 
phases. In the overlapping procedure of UV 1455A model 
light curves, we adopt the origin of time which is different 
from the origin of time in me of Table [3] This procedure is 
necessary for the same reason as in Figure [3] As a result, the 
peak brightness of each light curve is almost the same in the 
logt'-m' v plane as shown in Figure[6] In this figure, we adopt 
M env . = 1.5 x 10- 5 M Q for 1.05 M Q WD model (f s « 1.0) as 
a starting point of the model light curves, corresponding to 
m ff = 4. 1 of 1 .05 M Q WD model in Tabled 

We also obtain the absolute magnitudes of free-free model 
light curves in the same way as the CO novae mentioned 
in the previous subsections, but for the neon nova VI 974 
Cyg a s a comparison. A dopting the distance of d = 1.8 
kpc (IChochol et al. 1 119971) and the absorption of Ay = 1.0 
(IChochol etal. 1 11993) for VI 974 Cyg, we obtain the distance 
modulus of 

(m-M)y = [51og( £ //10)+Ay] VI974 Cyg = 12.3, (16) 

for V1974 Cyg. Then, we specify the absolute magnitude at 
the bottom of our free-free model light curves as 

M w = m' vl +2.5 log/ 8 - 12.3, (17) 

for V1974 Cyg and tabulate / s , m' w , and M w in Tabled We 
plot M w for both CO and Ne novae in Figure [5] This clearly 
shows that the brightness M w is smoothly increasing from ~ 5 
mag to ~ mag with Mwd from 0.55 to 1 .2 M for "CO nova 
2" (from 0.7 to 1.3 M Q for "Ne nova 2"). 

Now, we can determine the absolute magnitude of free-free 
model light curves, My = (mff - 15.0) + M W from Table [3] for 
various WD masses as shown in Figure|7] 

We also check whether or not My(15) is almost common 
among various neon novae. We have obtained My(15) = 
-5.6 ± 0.3 for 0.70-1 .3 M Q WDs (13 light curves with equal 
weight for all WD mass models), being rou ghly consistent 
with the empirical estimates (see Section [3~4i >. For the statis- 
tical point of view, neon novae have typical masses of 0.9- 
1.15 Mq. If we take this rage of WD masses (equal weight 
for 6 light curves), the value becomes My(15) = -5.85 ± 0.09 
as indicated in Figure [7] 

If we take all the CO and Ne nova light curves with equal 
weight (14+ 13 WD mass models), the average value becomes 
My(15) = -5.7 ±0.3. Thus, the absolute magnitudes of our 
free-free emission model light curves are consistent with the 
My(15) empirical relations. 

3.6. Consistency check among nova light curves 

We have proposed new methods for obtaining the absolute 
magnitude of novae in the previous subsections. Using these 
absolute magnitudes, we can estimate the distance to novae. 
In this subsection, we will check the accuracy of our new 
methods by comparing three nova light curves. 

First of all, we explain how to select a model that repro- 
duces well observational light curves. Figure [8] shows our 
method in the case of V1668 Cyg. We shift the model light 
curve "back and forth" and "up and down" to fit with the ob- 
servational y magnitudes. Together with the UV 1455A light 
curve, we finally determine by eyes the 0.95 M WD as the 



best-reproducing one among the light curves tabulated in Ta- 
ble |2] (0.55- 1.2 M by 0.05 M Q step). In this figure, we 
added two other light curves of different WD masses, i.e., 
0.85 and 1.05 M WDs. Here, we use different phases of the 
model light curves in Figure Q] to fit with the y observation; 
the 0.85 Mq WD model is shifted rightward (back) but the 
1.05 Mq WD model is shifted leftward (forth). As a result, 
the three model light curves of 0.85, 0.95, and 1 .05 Mq WDs 
give similar decline rates from maximum until day ~ 40. In 
other words, we can find similar shapes in different phases of 
different WD mass light curves in Figure [T] In this sense, it 
is hard to choose the best-reproducing model only from early 
phase optical and near IR data if neither UV 1455 A nor X- 
ray data are available. 

It should be also noted here that the ambiguity of WD mass 
estimate of + 0. 1 Mq introduces an error of about ± 0.4 mag 
in the distance modulus of (ra-M)y, when only early optical 
(or IR) light curves are available. In Figure[8] we obtain (m- 
M) v = m w -M w = 15.8-1.6 = 14.3 for Mwd = 0.95 Mq, but 
(m-M)y= 14.9-1.0= 13.9 for M WD = 1.05 M or (m-M)y = 
17.1-2.4= 14.7 for M WD = 0.85 Mq. 

Now let us go to the first example of a very fast nova 
V1500 Cyg (Nova Cygni 1975). Figure [9] demonstrates that 
optical V (filled triangles), y (filled circles), near infrared H 
(open squares), and K (open triangles) magnitudes are al- 
most overlapped with each other from ~ 5 days until ^100 
days after the outburst. This wavelength-free decline is a 
characteristic property of free-free emission. This impressive 
part is reproduced well by Equation (Q~|) with the models of 
Mwd = 1-05 M ("CO nova 2") and M WD = 1.15 M Q ("Ne 
nova 2"). V1500 Cyg is a superbright nova, in which the peak 
luminosity was about 4 mag bright er than t he Eddington limit 
of a 1 .0 Mq WD (e.g JFerland et al. Il986h . and the spectra of 
the superbright part ( first 5 days) are approxi mated by black- 
body emission (e.g., Gallaghe r & Nev 1119761) . Therefore, our 
free-free emission model cannot be applicable to the super- 
bright part. 

After day ~ 70, the V light curve gradually deviates from 
the other light curves (y, H, and K) because emission lines 
dominate the optical fluxes. The V magnitudes are 2 mag or 
more brighter than the y, H, and K magnitudes in the late 
phase . As extensively dis c ussed in Hachis u~& Kato I d2006l 
120071) and lHachisu et all (12008b , V magnitudes are often 
heavily contaminated by strong emi ssion lines such as [O III] 
A A 4959, 5007 (see lRead et aU2008l for V598 Pup spectra), 
which eventually dominate over the continuum, causing an 
increasing deviation from the free-free emission model. 

Comparing the apparent magnitude of VI 500 Cyg and the 
absolute magnitude of the model light curve of 1 .05 Mq (Ta- 
ble|2]), we can estimate the distance modulus. Figure|9]shows 
that the 1.05 M Q WD model is a best-reproducing model 
(magenta dash-dotted line) among other WD mass models 
in Table [2] for chemical composition of "CO nova 2." Using 
M w = 1.0 in Tableland m w = 13.0+0.5 = 13.5 from FigureHl 
we obtain the distance modulus to V 1500 Cyg, i.e., (m-M)y = 
m^-Mv, = 13.5- 1.0 = 12.5. If we adopt another chemical 
composition of "Ne nova 2," we obtain a best-reproducing 
model of 1.15 M Q WD (black solid line). Then we obtain the 
distance modulus of (m-M)y = m w -M w = 13.6-0.9 = 12.7 
with m w = 13.1 + 0.5 = 13.6 (Figure and M w = 0.9 (Table 
[3]). The difference between CO and neon nova models is only 
0.2 mag. Therefore, this method for obtaining the distance 
modulus is rather robust against the ambiguity of chemical 
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FIG. 8. — Close-up view of three light curves in Figures[T]and[2] Three models are for 1.15 (black thin line), 0.95 (blue thick line), and 0.85 Mq (Green thick 
line) WDs. None of them are squeezed or stretched. The model of 0.95 Mq WD reasonably reproduces both the optical y (and V) and UV 1455 A light curves of 
V1668 Cyg. In order to fit free-free model light curves of 1.15 and 0.85 Mq WDs with the observation, we shift the 1.05 Mq up and leftward, and the 0.85 Mq 
down and rightward. As a result, two free-free model light curves reasonably follow the y magnitudes but the two UV light curves do not fit with the observation 
of V1668 Cyg. Three arrows indicate to, ti, and (3 times of V1668 Cyg. The times of ?2 and (3 are determined along our model light curve of the 0.95 Mq WD, 
which we call "intrinsic" times. Our estimated "intrinsic" values of V1668 Cyg are ti = 14.4 and tj = 26 days. We also add observational V magnitudes (red 
small filled circles), which deviate from y magnitudes in the later phase due mainly to contribution from strong emission lines such as [O III]. Fortunately this 
deviation has no effects on the estimate of fj and (3 times because it starts about 40 days after the outburst in this case. Points A, B, and C correspond to the same 
initial envelope masses as in Figure |2]along the 0.95 Mq WD model. The magnitude at the end of a wind phase (outside the figure) is m w = 14.9, 15.9, and 17.1 
in this plot, for Mwd = 1 05, 0.95, and 0.85 Mq, respectively. So the distance moduli of these three light curves are (m—M)v = m w — M w = 14.9 — 1.0 = 13.9, 
15.8-1.6= 14.3, 17.1-2.4= 14.7, respectively. 



composition of a nova envelope. 

The distance to V1500 Cyg was estimated to be d = 1.5 ± 
0.2 kpc from the nebular (ejecta) expansion parallax (e.g., 
iSlavin et ai~lll995l). The ab s orptio n was obtained to be Ay = 
1 .60 ± 0.16bv lLance et al. I {[988), so that the distance mod- 
ulus to V1500 Cyg is (m-M) v = 12.5 ±0.4, being consistent 
with the value of 12.5 ("CO nova 2") and 12.7 ("Ne nova 2") 
estimated above based on our new method. 

Thus, the absolute magnitude estimates of VI 500 Cyg 
based on our free-free emission model light curves are con- 
sistent with the d istance derived fro m the nebular expansion 
parallax method (ISlavin et al. 1 1 19951) within an accuracy of 
0.2 mag in the distance modulus of (m-M)y. This in turn sup- 
ports our assumptions and simplifications on free-free emis- 
sion of novae. 

The second example is GK Per (Nova Persei 1901). We 
assume that, when two novae have a similar decline rate (or 
a similar timescale) of free-free light curve, their brightnesses 
are the same during the period in which the two free-free light 
curves are overlapped. Figure [10] shows such an example, a 
comparison between the light curves of GK Per and VI 500 
Cyg. GK Per has a decline timescale very similar to that of 
V1500 Cyg. Even in a typical transition phase (from day ~ 25 
until day ~ 150), we are able to nicely fit our universal decline 
law with the observation, especially the bottom line of each 
oscillation until day ~ 100 as shown in the figure. The top 
line connecting maxima of each oscillation deviates largely 
from our free-free model light curve. This deviation is caused 
mainly by the increase in continuum flux itself. The effect of 
strong emission lines starts from day ~ 100, where the bottom 



line of each oscillation starts to deviate from our model light 
curve. 

Our best-reproducing model is a 1 . 15 Mq WD among other 
WD masses in Table [3] Then we obtain the distance modulus 
of 

[(to-M) v ] ff = ot w -M w = 10.2-0.9 = 9.3 (18) 

from m w = 10.2 in Figure[10]and M w = 0.9 in Table[3] where 
FF means the method of our free-free (FF) light curve fitting. 

Direct comparison of the brightness between GK Per and 
VI 500 Cyg gives the difference of their distance moduli. 
From Figure \W[ we have 



[5lQg(<//10)+A v " 



VI 500 Cyg 



[51og(<//10)+Ay] GKPer = +3.4, 

(19) 

where d is the distance in units of pc and Ay is the absorption 
in the V band. Wit h d = 1 .5 kpc :and A v = 1 .6 for V1500 Cyg 
and d = 0.455 k pc (ISlavin et al. 119951) and A v = 3AE(B-V) = 
3. 1 x 0.3 = 0.9 dWuetal. II 198% for GK Per, we obtain 



(m-M) v = [51og(rf/10)+A v 



V1500 Cyg 



and 



(m-M) v = [51og(af/10)+A v 



I GK Per 



12.5, 



= 9.2. 



(20) 



(21) 



Therefore, the difference of 12.5-9.2 = 3.3 is consistent with 
the difference in Equation ( fT9b within an accuracy of 0.1 mag 
and the distance modulus of Equation d2Tb is also consistent 
with the result of the above FF method (9.3 mag) within 0.1 
mag. Thus, we confirm that the brightnesses of the two no- 
vae are the same within an uncertainty of ±0.2 mag. This 
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FIG. 9. — Light curves for VI 500 Cyg as an example of fbreak- Optical V 
(filled triangles), y (filled circles), near infrared (IR) H (open squares), and 
K (open triangles) magnitudes of the fast nova VI 500 Cyg. Light curves 
of free-free emission are essentially independent of wavelength, so that V, 
y, H, and K light curves of VI 500 Cyg are almost overlapped with each 
other from several days to ~ 70 days after the outburst. Then, the V light 
curve deviates from the others (y, H, and K) mainly because the wide V 
band is heavily contaminated by strong emission lines such as [O III]. Two 
other nova light curves of V598 Pup and V382 Vel show an almost identical 
decline from several days to several tens of days after the outburst, that is, 
the free-free part of light curves. O ptica l V (filled squares) and visual (small 
open circles) of V598 Pup (see Fig. ll6l for mor e de tails of V598 Pup) and V 
(open circles with plus) of V382 Vel (see Fig. 1 1 7 1 for more details of V382 
Vel). Black solid line: free-free model light curve of 1.15 Mq WD with 
a composition of "Ne nova 2." Magenta dash-dotted line: free-free model 
light curve of 1.05 Mq WD with a composition of "CO nova 2." Large 
open circle: the epoch when the wind stops. Dashed line: model light curve 
based on free-free emission following a law of F\ oc t~ 3 (Equation after 
winds stop). Near IR dat a of VI 500 Cyg are t aken from Ennis et al. ] 119771) . 
IGallagher &"Nev1 (WR$), and Kawar aetaTI (19761); y and V d ata of V1500 
Cyg are from Lockwood & Millis 1 1976) and Tempesti 11979). 
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FIG. 10. — Same as Figure [5] but we add the V light curve of GK Per 
(Nova Persei 1901) for comparison. We fit the smooth light curve of V1500 
Cyg with the bottom line of the GK Per light curve during the transition 
oscillations. Optical V data of GK Per are the same as those in Figure 2 of 
IHachisu & Katol BMfo . 

good agreement supports our assumption that two novae have 
the same brightness during their overlapping period if the two 
novae have a similar timescale as shown in Figure [7] 
The third example is V1493 Aql (Nova Aquilae 1999 



No.l). This nova shows a peculiar secondary peak (Fig- 
ure ITTb, which is seen in all t he bands, i.e., V, R, and I 
bands. Hachis u~& Kato I d2009l) explained this kind of sec- 
ondary peak by an additional energy source from strong mag- 
netic field (and rotation). Therefore, we regard that the nova 
light curves come back to a free-free emission light curve af- 
ter the additional energy input stopped. We see that, just after 
the secondary peak, the light curve sharply drops to merge 
into the light curve of V1500 Cyg. So we expect that V1493 
Aql has a similar decline rate to V1500 Cyg, i.e., our free- 
free model light curve of Mwd =1.15 Mq with "Ne nova 2" 
agrees well with the observation. The effect of strong emis- 
sion lines starts only from day ~ 100, where the V and visual 
magnitudes begin to depart from the / magnitude. 

Using a 1.15 M Q WD model in Table [3] we obtain the dis- 
tance modulus of 

[(m-M)y] FF = m w -M w = 19.3-0.9= 18.4 (22) 

from m w = 19.3 in FigureQUand M w = 0.9 in Tabled 

Direct comparison of the brightness between V1493 Aql 
and VI 500 Cyg gives the difference of their distance moduli 
as shown in Figure [TT] i.e., 

[51og( £ //10) + A v ] vl500 Cyg -[51og(,i/10) + A v ] VI493 Aql = -5.' 

(23) 

With (m-M) v = 12.5 from Equation (|20j for V1500 Cyg, we 
obtain the distance modulus of 

[(m-M)vkc = [51og(<//10)+A v ] VI493Aql = 18.2, (24) 

where LC means the method of our direct light curve (LC) 
fitting. The distance modulus of this nova was estimated to 
be (m-M)y = 17.8 dBonifacio et al. Il2 000b (m-M)v = 18.8 
dVenturini et al. 1 12004) , (m- M) v = 18.4 {A rkhipov a et al. I 
|2002D, and {m-M) v = 18.6 (iBurlakl [2008). based on the 
MMRD relations. The simple average values of these four 
observational estimates becomes (m-M)v = 18. 4 ±0.4, being 
consistent with our estimate of the FF and L C methods. If we 
adopt A v = 3AE(B-V) = 3.1 x 1.5 = 4.7 (lArkhipova et al. I 
120021) . the distance is d = 5.5 and 5.0 kpc for FF and LC 
method, respect ively. These values are c onsistent with d = 5.4 
kpc obtained by Arkhipova et al. (2002) based on the MMRD 
relations. 

In the present paper, we proposed two new methods for ob- 
taining distance modulus of a nova: one is fitting our free-free 
(FF) emission model light curve with the observation (we call 
this FF method). The other is fitting the light curve (LC) of 
a distance-unknown nova with the other distance- well-known 
nova (we call this LC method). From the above confirma- 
tions, we may conclude that our new methods give a reason- 
able distance. It should be addressed that our new methods for 
obtaining distance moduli of novae are applicable even when 
we missed the optical maxima of nova light curves. In Section 
[6] we apply these two methods to ten novae. 

4. CHARACTERISTIC TIMESCALES OF NOVA LIGHT CURVES 

There are several timescales that characterize nova light 
curves; to is the time in units of day at optical maximum 
(usually V magnitude); ti and ?3 are the times in units of 
day, in which V magnitude decays by 2 and 3 mag, respec- 
tively, from its maximum (from fo); fb rea k is the epoch when 
our free-free model light curves (or optical y, near infrared 
/, /, H, K light curves) bend sharply as shown in Figures [3] 
[6] and |9j f w i„d is the epoch when optically thick winds stop; 
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FIG. 1 1 .— Same as Figure|9] but we add light curves of V1493 Aql 1999#1 
for comparison. We fit the smooth light curve of VI 500 Cyg with the bottom 
line of the V1493 Aql light curve. Optical V (filled squares), R (open circles 
with a plus), I (open diamonds), and visual (very small open circles) are taken 
from the American Association of Variable Star Observers (AAVSO) and the 
Variable Star Observers League in Japan (VSOLJ). Some visual magnitude 
estimates are also taken from IAU Circulars 7225, 7228, 7232, 7258, 7273, 
and 7313 (open squares). Here we assume the outburst day (day 0) of foB = 
JD 2451369.0 (1999 July 9.5 UT) for V1493 Aql. 

and fH-hnrniny is the epoch whe n hy drogen shell-burning e nds. 
Hachisu&Kato (20061 12007b and lHachisu"eTaI1(l2008h cal- 
culated free-free model light curves for various WD masses 
with many sets of chemical composition. They fitted their 
model light curves with observation and derived fbreak, fwind, 
and fH-buming for many novae. Table |4] summarizes such val- 
ues of fbreak, fwind, and ^-burning f° r 15 classical novae, taken 
from our published papers as well as £2 and £3 for the same 15 
novae, which are taken from literature. Figures [T2T4T41 show 
these five timescales for the 15 novae. 

When the optically thick wind stops, the photospheric ra- 
dius shrinks drastically and th e photospheric tempera ture in- 
creases up to 20-30 eV (e.g. jHachisu & Ka to 2003) so that 
the supersoft X-ray flux emerges (see Appendix B for the pho- 
tospheric radius which is shrinking with the temperature of 
the WD envelope during a supersoft X-ray phase). We regard 
that f w i n d is almost equal to the emerging time of supersoft 
X-ray, fx-on- After the hydrogen shell-burning stops, the WD 
envelope quickly cools down and the supersoft X-ray phase 
ends. So we regard that fH-buming is almost equal to the decay 
time of supersoft X-ray, fx-off- Among these five timescales 
of ?2, £3, fbreak, fwind, and fH-buming, the first three timescales, 
f2, f3, and fbreak, are much shorter than the other two, f w ; n( i and 

f H— burning ■ 

4.1. Relations between various timescales 

In what follows we introduce several relations between 
these five timescales, f 2 , £3, fbreak, fwind, and fH-buming- In Fig- 
ure H2 we plot fH-buming, fwind, and f bre ak against fbreak, which 
are the oretical ones calculated f rom o ur free-free model light 
curves (Hachis u & Kato Il2006l 12007b . These values depend 
on Mwd and the chemical composition, but are actually in- 
dependent of M env 0- This is because the difference in M e nv,o 
makes a large difference in the optical maximum but affects 
only the very early phase of nova outbursts as shown in Fig- 
ures [2] and [8] and the difference in the epoch of optical max- 
imum is so short (a few to several days) that it can be ne- 
glected compared with the timescales of fH-buming, fwind, and 



fbreak (several tens of days or a few hundred days). So we 
can derive relations among fH-buming, fwind, and fb rea k from our 
models (independently of observations) as shown in Figures 
[T214T41 Once we determine fbreak of individual novae from 
observations, we can theoretically predict the duration of a 
luminous supersoft X-ray phase from Figure [T21 
From Figure [T2l we obtain approximate relations between 

fbreak, tUm-On (f X -on ~ fwind), and tumoff (fx-off ~ fH-buming) 

times as 

fx-on ~ fwind = (2-6 ± 0.4)f b reak, (25) 



fx-off ~ fH-buming = (0.71 ± 0.1) (fbreak) 



and 



1.5 



fx-off « fH-burning = (0. 17 ± 0.023)(f wind ) 1 5 , (27) 



(26) 



for 25 < fbreak ^ 250 days. 

It should be noted that the proportionality of f w ; n d to /break, 
i.e., f„ind/fbreak ~ 2.6, is a direct result of our scaling law (the 
universal decline law), because both fb rea k and f wrn d are during 
the optically thick wind phase which is the subject of univer- 
sal decline law. Then the f w i nr j is directly proportional to fbreak- 
On the other hand, the turnoff time of hydrogen shell burning, 
fn-buming, is not proportional to fbreak or f w ind, because the dura- 
tion of static hydrogen shell burning phase depends on hydro- 
gen content X (amount of fuel) of the envelope and WD mass 
(temperature) in a different way from the optically thick wind 
phase, in which the decay rate is determined mainly by wind 
mass-loss not by hydrogen burning rate (see Figure[6]). There- 
fore, Equations d26b and d27i > are empirical laws obtained nu- 
merically, the power of which is not 1.0 but approximately 
1.5. 

It is, however, not easy to obtain fb rea k because nova light 
curves such as V magnitudes are often heavily contaminated 
by strong emission lines in the later phase and this contami- 
nation easily clouds the break point of light curves as can be 
seen in Figure [9] In the next subsection, we derive relations 
between £3 time and two of f w j n d and fH-buming times, because 
£3 time is determined relatively easily from early optical ob- 
servation. 

4.2. Prediction formulae by intrinsic t?, time 

It is very useful if we can predict fx- on (~ fwind) and fx-off 
(ps fH-buming) from £3 time, because £3 is popular and widely 
used to define the nova speed class. In this subsection, we 
derive relations between £3 and fx-on, and between £3 and fx-off- 

As already explained in the previous subsection, fbreak, fwind, 
and fH-buming are practically independent of the initial enve- 
lope mass, M env .o, but £3 time depends on the peak brightness 
of an individual nova. As shown in Figures |2] and [HJ the ini- 
tial envelope mass M env .o is larger, the optical maximum is 
brighter. Therefore, £2 and f3 depend on M env .o- This is the 
reason why we cannot theoretically derive f 2 and f3. We need 
information of maximum brightness to specify £3 time for in- 
dividual novae. 

Now we go back to the well-observed VI 668 Cyg which 
we regard as a typical classical nova. Thus, we adopt point B 
in Figures|2]and[8]as a typical (an average) case of maximum 
brightness. In this case, £2 = 14.4 and £3 = 26 days are de- 
rived from Figures|2]and[8] Then, we have f3/fbreak = 26/100 
(see Table|6]for V1668 Cyg timescales). Using the universal 
decline law, we obtain a relation between £3 and fb re ak as 



£3 



h/fs 



26 ±4 



fbreak fbreak / fs 100 



= 0.26 ±0.04, 



(28) 



Hachisu & Kato 



TABLE 4 

Light Curve Parameters of Novae in Literature 



obj cct 




'3 


ret. a 


'break 


'wind 


'H— burning 


ref. b 






tdavs) 






tdavs) 


( davs) 




GKPer 1901 


7 


13 


1 


73 


182 


382 


5 


V1500Cyg 1975 


2.9 


3.6 


2 


70 


180 


380 


8 


VI 668 Cyg 1978 


12.2 


24.3 


3 


110 


280 


720 


8 


V1974 Cyg 1992 


16 


42 


1 


96 


250 


600 


8 


V382 Vel 1999 


6 


12 


4 


49 


120 


220 


9 


V2361 Cyg 2005 


6 


8 


5 


69 


169 


340 


5 


V382 Nor 2005 


12 


18 


5 


73 


182 


382 


5 


V5115 Sgr2005 


7 


14 


5 


60 


145 


280 


5 


V378 Ser 2005 


44 


90 


5 


257 


858 


2560 


5 


V5116Sgr2005#2 


20 


33 


5 


114 


319 


757 


5 


V1188Sco 2005 


7 


21 


5 


44 


110 


190 


5 


V1047 Cen 2005 


6 


26 


5 


257 


858 


2560 


5 


V476 Set 2005 


15 


28 


6 


100 


260 


590 


5 


VI 663 Aql 2005 


13 


26 


5 


100 


260 


590 


5 


V477 Set 2005#2 


3 


6 


7 


32 


80 


121 


5 



3 tz and t$ times at observational face value taken from literature of l{Downes & Duerbeck 
(2000). 2-Warner (1995), 3 Mallama & Skillman (1979), 4- Vanlandinsham et al. (2001), 5- 
Hachisu & Kato i 2007), 6 Munari et al. i 2006b), 7 Munari et al. :2006a) 

b 'break, 'wind, and 'H-buming taken from 5£Hachisu & Kato i 2007), 8 Hachisu & Kato (200§), 9- 
present work 



TABLE 5 

Prediction Formulae vs. Nova Decay Timescales in Literature 



object 


h 

. (day) 


(day) 


t 2 -t 3 
(eq.(29]) 


'3 -'Xon 

(eq.HD) 


'3 ~'Xoff 
(eq.CED) 


reference" 


V598 Pup 2007#2 














V382 Vel 1999 


. 7.5 


15 


<yes> b 




no 


Burlak (2008) 


V382 Vel 1999 


. 4.5 


9 


<yes> 




<yes> 


Delia Valle et al. (2002) 


V382 Vel 1999 


. 6 


12 


yes 




yes 


Vanlandingham et al. (200 1 ) 


V382 Vel 1999 




12.5 






yes 


Liller& Jones (2000) 


V382 Vel 1999 


. 6 


10 


yes 




yes 


Delia Valle et al. (1999) 


V4743 Sgr 2002#3 . 


. 12 


22 


yes 




<yes> 


Burlak (2008) 


V4743 Sgr 2002#3 . 




15 






<yes> 


Nielbock & Schmidtobreick (2003) 


V4743 Sgr 2002#3 . 


. 9 


16 


yes 




yes 


Morgan et al. (2003) 


VI 281 Sco2007#2 . 














V597 Pup 2007#1 


'. 2.5 










Naiketal. (2009) 


V1494 Aql 1999#2 . 


. 8 


16 


<yes> 


no 


no 


Burlak, (2008) 


V1494 Aql 1999#2 . 


. 7 


23 


no 


yes 


<yes> 


Venturini et al (2004) 


VI 494 Aql 1999#2 . 


. 6.1 


15.8 


no 


no 


no 


Barsukova & Goranskii (2003) 


VI 494 Aql 1999#2 . 


. 6 


15 


no 


no 


no 


Arkhipova et al. (2002) 


V1494 Aql 1999#2 . 


. 6.6 


16 


no 


no 


no 


Kiss & Thomson (2000) 


V2467 Cyg 2007 


. 7.3 


15.1 


<yes> 




no 


Lynch et al (2009) 


V2467 Cyg 2007 


. 7.6 


14.6 


yes 




no 


Poggiarii^009) 


V2467 Cyg 2007 


. 8 










Senziani et al. (2008) 


V5116Sgr2005#2 . 


. 6.5 


20.2 


no 




no 


Dobrotka et al. (2008) 


V5116Sgr2005#2 . 


. 18 


45 


no 




no 


Burlak~(2Tj0c!r 


V5116Sgr2005#2 . 


. 20 


33 


yes 




yes 


Hachisu & Kato (2007) 


V574 Pup 2004 


. 37 


85 


no 


no 


no 


Burlak (2008) 


V574 Pup 2004 


. 13 


58 


no 


no 


no 


Siviero etaT (2005) 


V458 Vul 2007 


. 8 


31 


no 


no 




Wesson etal. (2008) 


V458 Vul 2007 


. 7 


15 


<yes> 


no 




Poggiani (2008) 


V458 Vul 2007 


. 7 


18 


no 


no 




Tarasova (2008) 


VI 974 Cyg 1992 


. 22 


48 


no 


no 


no 


Burlak (2008) 


V1974 Cyg 1992 


. 16 


42 


no 


no 


no 


Warner (1995) 


V1668 Cyg 1978 


. 12.2 


24.3 


yes 






Mallama & Skillman (1979) 


V1500 Cyg 1975 


. 2.9 


3.6 








Warner (1995) 



reference for quoted values of tj and at observational face value 

"yes" satisfies eqs.{29), D0K and \3\ \ but "no" does not and parenthesis means slightly outside the border 
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FIG. 12. — Five characteristic nova timescales against (break- From top 
to bottom, hydrogen shell-burning (/H-buming). optically thick wind (f wln d), 
break point of nova light curve (fbreak). '3. and h- Supersoft X-ray on/off 
correspond to ? w i n d and fH-bumine, respectively. Each symbol in the upper left 
corner represents o ur various WD models with different chemical composi- 
tions taken from Hachisu & Kato 12006). We also plot observational face 
values of tj {filled triangles) and ?2 (open circ les) for 15 novae in Table [4] 
Thick red solid line: central values of Equati on 128) . Upper/lower red dashed 
lines: upper/lower limit values of Equation )28t . Thick blue solid line: cen- 
tral values of f2-fbreak relation. Upper/lower blue dashed lines: upper/lower 
values of ?2 -fbreak relation. Colored marks indicate t% and (3 times that devi- 
ate largely from the statistical relations. Green: visual magnitudes of VI 188 
Sco are so contaminated by emission lines that is largely deviated from the 
central line. Red: visual and V light curves of V2361 Cyg are affected by 
dust formation. Blue: V1500 Cyg is a superbright nova and its face values of 
fj and f3 times do not follow the universal decline law. Magenta: visual and 
V light curves of VI 047 Cen is also a superbright nova. 

for 8 < f3 < 80 days, where the upper and lower bounds 
are corresponding to ±0.5 mag brighter/darker (more mas- 
sive/less massive initial envelope mass) points than that of 
point B. This relation is indicated by a red thick solid line in 
Figure [12] together with the upper and lower red dashed lines 
that denote upper/lower values (±0.04) of Equation d28l ). 

There is a relation between f 2 and £3 derived from the scal- 
ing law of o ur universal decline law, which is also well used 
in literature. Hachisu & Kato] d2007l) have already derived 



r 2 = (0.6±0.08)r 3 , 



(29) 



for the optical light curves that follow the free-free model light 
curve with 8 5, h < 80 days. This re lation has been empiri- 
cally known (Capaccioli et al. 1 1 19901) and theoretically con- 
firmed by Hachisu & Kato (2006) based on the universal de- 
cline law with a slope of F\ oc f" 1 75 as shown in Figures[3]and 
[6] Inserting the central value of Equation (|29l i into Equation 
d28l ). we obtain relations between fb rea k and ti relations by blue 
solid and dashed lines as shown in Figure [121 
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FIG. 13. — Same as Fi gure [121 but against t w - m( i. Thick red solid line: 
central values of Equation )30t . Uppe r/lower red dashed lines indicate up- 
per/lower limit values of Equation )30t . Thick blue solid line: central values 
of *2-f W md relation. Upper/lower blue dashed lines: upper/lower values of 
f 2 _f wmd relation. See text for more details. 

Figure Q~2] also shows f 2 and times for individual novae 
in Table [4] Four novae indicated by colors are far from the 
relation of Equation ( |28l , but there is a special reason in each 
nova as described in the figure caption. The other novae rea- 
sonably follow the above relation with the reasonable scatter 
corresponding to ±0.5 mag brighter/darker of the peak bright- 
ness. This supports that V1668 Cyg has an average M env o. 

In the same way, we derive a relation between ?3 and f w i nt j 
from Equations d25l > and d28l , i.e., 



tx- 



^wind = (10.0±1.8)f 3 , 



(30) 



for 8 < ?3 < 80 days. Figure Qj] shows the five timescales 
against f w ; n( i and we also see a reasonable fit of our relation 
(Equation d30l l) with the observation. 

From Equations d26l i and d28l . we derive a relation between 
f3 and fH-buming, i.e., 



f X-off ~ fH-burning - (5.3 ± 1 .4)(f 3 ) 



1.5 



(31) 



for 8 < ?3 < 80 days. Figure [14] shows the five timescales 
against fH-buming and we can also see a reasonable fit of Equa- 
tion OTb with the observation. 

However, it is not always easy to determine % (or t£) time 
from observation. In fact, Table [5] shows nova timescales ti 
and f3 of 1 1 novae taken from literature, which are sometimes 
very different from author to author, and sometimes very dif- 
ferent from our estimates in Table |6] Therefore, we propose 
a more robust way. Instead of observational face values of ti 
and ?3, we will use "intrinsic" ti and ?3, which are obtained 
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FIG. 14. — Same as Figu rell2l but against ^H-buming- Thick red solid line: 
central values of Equation J3U . Uppe r/lower red dashed lines indicate up- 
per/lower limit values of Equation )3U . Thick blue solid line: central values 
of /2-'H-burning relation. Upper/lower blue dashed lines: upper/lower values 
of ?2-'H-burning relation. See text for more details. 

from our model light curves fitted with observation. For ex- 
ample, we use ?2 = 14.4 and f 3 = 26 days for V1668 Cyg ob- 
tained from Figures |2] and [8] instead of observational face val- 
ues. This "intrinsic" £2 and £3 times are especially useful when 
the observational data deviate from our model light curve near 
the peak. In what follows, we will see how to determine "in- 
trinsic" ?3. 

Figure[8]explained how to select our best model light curve 
from early phase observations. It also shows that three dif- 
ferent WD mass (0.85, 0.95, and 1 .05 M Q ) models can repro- 
duce the early phase light curve and we can hardly select one 
only from the early phase optical and IR light curves. How- 
ever, this uncertainty of mass determination is not a problem 
in the estimation of "intrinsic" f 3 time, because these three 
model light curves, of course, give similar values of £ 3 time 
(f 3 ~ 26 ± 2 days) along our model light curves. 

Using relations of Equations (|2"5ll-(l3Tll, we can predict the 
period of a supersoft X-ray phase of individual novae from £3 
(or ?2, fbreak) time. In Section|6] we will show that our formu- 
lae are much better fitted with observation if we use intrinsic 
£ 3 time than if we use face values of £3. Therefore, in what 
follows in this paper, we use intrinsic £3 time in order to pre- 
dict f X -on ~ fwind and fx-off ~ fH-burning from Equations ([30]l 
and CD}. 

We assumed here that a nova has a typical M env 0, i.e., has 
a peak brightness similar to that of V1668 Cyg (i.e., point B 
in the rescaled t/f s -m v plane, Figure|2]i. If the rescaled peak 
magnitude is much brighter/darker than that of VI 668 Cyg 



(for example, points A/C in Figure|2]l, X-ray on and off times, 
fx-on w fwind and fx-off « fH-burning, can deviate from the predic- 
tion of Equations ( f30b and OTb . In this sense, our prediction 
formulae (f30b and OTb are statistical ones. 

5. MMRD RELATION DERIVED FROM A UNIVERSAL DECLINE 

LAW 

The Maximum Magnitude-Rate of Decline (MMRD) re- 
lation i s often used to estimate the distan c e to novae (e.g., 
I Cohen l fl98l iDella Valle & Livio I H995t I de Vaucouleurs I 
Il978t iDownes & Duerbeck 1 120001; ISchmidtlll957l) . In this 
section, we derive theoretical MMRD relations based on our 
universal decline law. 

5.1. Average MMRD relation 

Figures |4] and [7] clearly shows a trend that a more massive 
WD has a brighter maximum magnitude (smaller My max ) and 
a faster decline rate (smaller £3 time). The relation between 
f3 and Mv,max for novae is called "the Maximum Magnitude- 
Rate of Decline" (MMRD) relation. In order to derive theo- 
retical MMRD relations, we go back to Figure |2] and utilize 
the universal decline law starting from point B (a typical case 
of nova envelope mass). 

We have already shown that all the free-free emission model 
light curves for different WD masses overlap with each other 
(they converge to 0.95 M Q WD in Figure [2] against V1668 
Cyg data) when they are squeezed in the direction of time by 
a factor of / s , data of which are tabulated in Table [2] Now 
we assume that all the free-free emission model light curves 
reaches the same maximum point, like point B in Figure [2] 
In such a case, the apparent maximum brightness my .max of 
each light curve with different WD masses is expressed as 
wJv im ax = OT !/max + 2.51og/ s . The tj, time of each model light 
curve with different WD masses is squeezed to be ? 3 = f s ty 
Eliminating / s from these two relations, we have 

m v ,ma = 2.51ogf3 + m( /rnax -2.51ogf3. (32) 

Since point B corresponds to the optical peak of V1668 Cyg 
{f s = 1 in Figure 0, we have already known ti = 26 days and 
m 'vmax = 6-2- Substituting these values of (fj = 26,m' Vmax = 
6.2) into Equation d32l . we have 

'«y,max = 2 .5 log f 3 + 6 .2 - 2 .5 log 26 = 2 .5 log 1 3 + 2 .7 . (33) 

This is our theoretical apparent MMRD ration for a starting 
point B in Figure [2] Next, using the distance modulus of (m— 
M)y, we obtain our theoretical MMRD relation as 

My.max = max ~ ('« ~M)y 

= 2.5 logf 3 + m' Vm . dx - 2.5 log t' 3 - {m-M) v . (34) 

Substituting (m— M)y = 14.3, the distance modulus of V1668 
Cyg, into Equation d34"i >. we have 

My max = 2.5 logf 3 + 6.2 - 2.5 log26 - 14.3 = 2.5 log f 3 - 11.6. 

(35) 

Figure [15] shows this theoretical MMRD relation for novae 
that has a typical (average) M env ,o like V1668 Cyg. In the right 
axis we convert the apparent magnitude to the absolute mag- 
nitude by using the distance modulus of V1668 Cyg. Point B 
in Figure [TBI corresponds to V1668 Cyg (f s = 1: 0.95 M Q , in 
which nova light curves start from point B in Figure |2])- The 
other WD mass cases starting from the same point B in Fig- 
ure |2] have different / s ^ 1. Changing / s (or M W dX we have 
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FIG. 15. — The maximum magnitude-rate of decline (MMRD) relation against t-j time based on the universal decline law. Point B represents a nova on a 
O.95M0 WD with a typical initial envelope mass. The black solid line passing through point B is corresponding to an MMRD relation of novae having a typical 
initial envelope mass but different WD mass. The black solid line passing through point A is an MMRD relation for a much larger initial envelope mass than that 
of point B (much lower mass accretion rates). The black solid line passing through point C is an MMRD relation for a much smaller initial envelope mass than 
that of point B (much higher mass accretion rates). Therefore, the region between upper and lower black solid lines represents our theoretical MMRD relation. 
Green solid lines connect the same WD mass, from 1.15 to 0.65 Mq by 0.1 Mq step (green attached number) for "CO nova 2," but from 1.25 to 0.75 Mq by 
0.1 Mq step (blue attached number) for "Ne nova 2." The absolute magnitude My, raax at maximum is obtained from calibration with VI 668 Cyg of the distance 
modulus (m-M)v = 14.3. Red filled star marks: the recurrent nova RS Oph as an example fo r a very high mas s accretion rate and very short timescale ?3 (very 
small f s ). Two well-known MMRD rela tion s are added: Kaler-Schmidt's law (Schmidt 1957) of Equation (38) labeled "MMRD1" (blue solid lines), and Delia 
Valle & Livio's (1995) law of Equation (39) labeled "MMRD2" (magenta solid lines). We add upper/lower bound s of ±0.5 m ag for each of these two MMRD 
relations. Red filled circles: We also add MMRD points for individual novae taken from Table 5 of Downes & Duerbeck (2000). 



a line passing through point B in Figure [15] This is our the- 
oretical MMRD relation, the parameter of which is Mwd (or 
/ s ). We also show the equivalent parameter, stretching factor 
f s (against point B), in the upper axis of the same figure. 

The black thick solid line passing through point B shows 
quite a good agreement with the central blue and magenta 
lines, which are empirical MMRD relations obtained from ob- 
servation as will be explained bellow. 



5.2. Scatter of MMRD relations from average 

If the initial envelope mass M env ,o is larger or smaller than 
that of point B (which is an average value, see Figure |2}, we 
obtain different MMRD relations. In Figure [15] we show five 
points of green filled squares corresponding to different initial 
envelope mass of M env ,o- Point A (t^ = 16 days, m' Vm!ix = 4.6) 
corresponds to a larger envelope mass at ignition (point A in 
Figure |2j. The black solid line passing through point A is 
other »?y, m ax-logf3 relation for a much larger initial envelope 



mass of Menv.o ( at point A in Figure[2]i, i.e., 

M Vmax = 2.5 log? 3 +4.6-2.51ogl6- 14.3 = 2.51ogf 3 - 12.7. 

(36) 

On the other hand, the black solid line passing through point 
C (fj = 43 days, m' Vmllx = 7.9) is also different m v . mdx -\og £3 
relation for a much smaller initial envelope mass (at point C 
in Figure|2]i, i.e., 

M 1/max = 2.5 logf 3 + 7.9-2.5 log43- 14.3 = 2.5 logf 3 - 10.5. 

(37) 

For a given WD mass, M env ,o is larger (smaller) for a smaller 
(larger) accretion rate. Points A, B, and C correspond to dif- 
ferent ignition masses for a given WD mass. When the accre- 
tion rate is much smaller than that for point B, the brightness 
reaches, for example, point A at maximum and then we have a 
shorter f 3 time (see Figure©. In the case of a larger accretion 
rate than that for point B, the brightness goes up only to point 
C at maximum and we have a longer £3 time. So, even for 
a given Mwd and chemical composition of the envelope, we 
have different f 3 time, depending on M env ,o- In other words, a 
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variation of mass accretion rates from an average value makes 
a divergence (or width) in the MMRD relation for individual 
novae. 

The five points including points A, B, and C are calculated 
for VI 668 Cyg model light curves, all of which correspond to 
the 0.95 Mq WD of "CO nova 2" (/ s = 1). The green solid 
line connecting these five points represents the My max -log t^ 
relation for 0.95 M Q WD of "CO nova 2." We then obtain 
other My max -logf3 lines for other WD masses as shown by 
green lines (attached numbers represent WD masses) in Fig- 
ure[15] In the case of neon novae, we also obtain similar lines, 
which are also shown in Figure [15] These lines are almost 
overlapped to those of "CO nova 2" but different WD masses 
(blue attached numbers for WD masses). 

5.3. Comparison with other MMRD relations 

Our MMRD (My max -logf3) relation passing through point 
B is in a good agreement with two well-known MMRD 
relations empiri cally obtained: one is Kaler- Schmidt's law 
(ISchmidtl [T957). i.e., 



My. max = -11.75 + 2.51ogf 3 , 
and the other is Delia Valle & Livio's (1995) law, i.e.. 

1.32-log^ 



Mi 



V.max 



-7.92-0.81 arctan 



0.23 



(38) 



(39) 



where My max is the absolute V magnitude at maximum and 
we use Equation d29l to calculate ?2 from in Equation ( 1391 
in Figure [15] Kaler-Schmidt's law is denoted by a blue solid 
line with two attendant blue solid lines corresponding to ±0.5 
mag brighter/darker cases. Delia Valle & Livio's law is indi- 
cated by magenta solid lines. These two well-known empiri- 
cal MMRD relations are very close to our theoretical MMRD 
relation at average. 

We also add observational MMRD points (My max and 
tj) for i ndividual novae, data of whi ch are taken from Ta- 
ble 5 of iDownes & Duerbeck I (120001) . although their t^ are 
not intrinsic but observational face values. It is clearly 
shown that a large scatter of individual points (red filled 
circles) from the two empirical MMRD relations falls into 
between upper/lower cases of our MMRD relations for the 
largest/smallest initial envelope mass of M env o- This simply 
means that there is a second parameter to specify the MMRD 
relation for individual novae. The main parameter is the WD 
mass represented by stretching factor / s . The second param- 
eter is the initial envelope mass (or the mass accretion rate to 
the WD). This second parameter can reasonably explain the 
scatter of individual novae from the empirical MMRD rela- 
tions ever proposed. 

It is interesting to place the position of the recurrent nova 
RS Oph in this diagram. RS Oph locates in the smallest cor- 
ner of ?3 time and a relatively darker level of My max , near 
the My max -logf3 relation passing through point C. This in- 
dicates a much smaller envelope mass at the optical maxi- 
mum, suggesting that the mass accretion rate to the WD is 
very high. This situation is very consistent with the total 
picture of recurrent novae; a very massive WD close to the 
Chandrase khar mass and a high mass accretion rate to the 
WD (e.g.. lHachisu & Kato 11200 ll) . In this figure, we adopt 
My. ^ = -8.3 and t-j = 10.5 days with the distance of d = 1.6 
kpc dHiellming et al. Ill986h absorption of Ay = 3.1E (B-V) = 
3Tx 0.73 = 2.3 dSniiders1ll987allbh . m Viinax = 5.0 dRosinol 
11987b . and t^ = 10.5 days from optical light curve fitting 
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FIG. 16. — Optical and supersoft X-ray light curves for V598 Pup. We plot 
free-free model light curves (labeled "Opt&IR") and 0.2 — 0.6 keV supersoft 
X-ray fluxes (labeled "X-ray") of three WD mass models of Mwd = 1-3 Mq 
(dash-dotted line), 1.28 Mq (thick solid line), and 1.25 Mq (thin solid line), 
for the envelope chemical composition of "Ne nova 2." We select a best 
model of 1.28 Mq WD among three WD mass models, in order to reproduce 
the supersoft X-ray data. The supersoft X-ray data is sandwiched between 
the 1 .3 and 1 .25 Mq WDs, so we safely regard that M WD = 1 .28 ± 0.03 Mq . 
The X-ray flux data (lar ge diamonds) without error bars are taken from 
Read et al. I f2 007. 2008). Optical V data of V598 Pup are taken from 
Pojmanski et al. (2007) and visual data are from the American Association 



of Variable Star Observers (AAVSO). Arrows indicate two epochs when the 
optically thick wind stops about 90 days after the outburst and when the 
steady hydrogen shell burning atop the WD ends about 150 days after the 
outburst for our 1.28 Mq WD model. We added optical and near IR data of 
V1500 Cyg for comparison (see Figure[9j- 

with our free-free model light curves (lHachisu & Kato ll200lt 
lHachisu et al. Il2006l I2007T) . 

6. LIGHT CURVE ANALYSIS OF TEN NOVAE WITH SUPERSOFT 
X-RAY EMISSION 

In this section, we examine ten classical novae which we 
have not yet analyzed in our series of papers, i.e., V598 Pup, 
V382 Vel, V4743 Sgr, V1281 Sco, V597 Pup, V1494 Aql, 
V2467 Cyg, V51 16 Sgr, V574 Pup, and V458 Vul, in the or- 
der of increasing timescale. All are supersoft X-ray detected 
novae, so we can check whether or not our prediction formu- 
lae of a supersoft X-ray phase work on these novae. 

6.1. V598 Pup 2007#2 

V598 Puppis (Nova Puppis 2007 No. 2) was ser endipitously 
disco ve red by the XMM-N ewton slew survey dRead et al. I 
2008). Tor res et al. I (|2007l) subsequently found an optically 
bright counter part of this X-ray object. The spectra obtained 
on November 16.34 UT (day ~ 170) reveal numerous features 
suggesti ve of a nova that already entered the auroral phase 
(see also lRead et al. 1120081) . Analysis of All Sky Automated 
Survey (ASAS) archive showed that the V magnitude rose to 
my = 4.1 from my > 14 between June 1.966 UT and June 
5.968 UT (Poimansk i et al. II2007I) . Since we cannot specify 
the exact date of the outburst, we assume here that the outburst 
day is ? 0B = JD 2454253.0 (2007 June 1.5 UT). 

The X-ray satellite Swift detected X-ray flux from V598 
Pup, which declined rapidly from Fx(0.2-2keV) = 1.5 x 
10" 10 (day ~ 130) to 6.7 x 10" 12 (day - 150) and then to 
1.3 x 10~ 12 erg s" 1 cm -2 (at day ~ 170 after the outburst) 
dRead et al. 1 120071) as shown in Figure [16] This quick drop 
in the X-ray flux clearly showed that hydrogen burning on the 
WD ended around 150 days after the outburst. 

Figure[T6]shows observed optical light curves of V598 Pup 
as well as the X-ray data. Unfortunately there are only two 
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TABLE 6 

TlMESCALES OF NOVAE ESTIMATED FROM UNIVERSAL DECLINE LAW 



object 


■ f2 a 
(day) 


?3 a 
(day) 


'break 
(day) 


fwind 
(day) 


*H— burning 
(day) 


^emerge 

(day) 


M W D 

(Mq) 


chem.comp. b 


VDVo rup 2(JU/#2 


~ , 






yo 


1 50 




1.28 


Ne nova 2 


Anon i f^i i Ann 

v3o2 vei iyyy 


6.6 


12.4 


52 


120 


220 


32 


1.23 


Ne nova 2 


V 474.5 Sgr 20(j2#3 


9.4 


17 


70 


180 


380 


36 


1.15 


Ne nova 2 


V 12ol oco 2UU/#2 


9.5 


1 £ "7 

16. / 


78 


1 90 


425 




1.13 


Ne nova 2 


V597 Pup 2007#1 


. 8.7 


16.5 


74 


200 


430 


52 


1.2 


Ne nova 3 


V1494 Aql 1999#2 


11.1 


21.2 


96 


260 


670 


66 


1.13 


Ne nova 3 


vi4y4Aqi iyyy#2 


11.7 


22.1 


102 


260 


670 


72 


1.06 


Ne nova 2 


v I4y4 Aqi lyyyffz 


1 A 1 

14. / 


ZO.D 


1 20 


285 


ID 


95 


o.yz 


CO nova 2 


V2467 Cyg 2007 


. 14.3 


25.3 


110 


290 


765 


68 


1.11 


Ne nova 3 


V2467 Cyg 2007 


14.8 


27.9 


113 


290 


765 


75 


1.04 


Ne nova 2 


V2467 Cyg 2007 


17.6 


31.7 


135 


310 


750 


96 


0.90 


CO nova 2 


V5116Sgr2005#2 


16.7 


30.4 


127 


355 


1000 


91 


1.07 


Ne nova 3 


V5116Sgr2005#2 


16.2 


30.5 


134 


350 


995 


100 


1.0 


Ne nova 2 


V5116Sgr2005#2 


20.0 


35.8 


169 


390 


1010 


130 


0.85 


CO nova 2 


V574 Pup 2004 


18.6 


35.0 


144 


390 


1140 




1.05 


Ne nova 3 


V458 Vul 2007 


20.4 


37.2 


145 


415 


1250 


62 


0.93 


CO nova 4 


V1974 Cyg 1992 


16 


26 


108 


250 


600 


90 


1.05 


Ne nova 2 


V1668 Cyg 1978 


14.4 


26 


100 


280 


720 


87 


0.95 


CO nova 3 


V1500 Cyg 1975 


. 7.2 


13 


70 


180 


380 


50 


1.15 


Ne nova 2 



a ti and ?3 times are all intrinsic values estimated along our theoretical light curve 
b see TableQ] 



TABLE 7 

Summary of Physical Properties of Novae 



object 


■ Mwd 


M 2 


Porb 


a 


R* 


M wmd a 


Menv.O b 


chem.comp. c 




(Mq) 


(Mq) 


(day) 


(«©) 


(«©) 


(1O- 5 M ) 


(1O- 5 M ) 




V598 Pup 2007#2 


. 1.28 










0.31 


0.34 


Ne nova 2 


V382 Vel 1999 


. 1.23 


0.31 


0.1461 


1.35 


0.68 


0.43 


0.48 


Ne nova 2 


V4743 Sgr 2002#3 . 


. 1.15 


0.70 


0.2799 


2.2 


0.93 


0.69 


0.78 


Ne nova 2 


V1281 Sco2007#2 . 


. 1.13 










0.92 


1.02 


Ne nova 2 


V597 Pup 2007#1 


. 1.2 


0.22 


0.1112 


1.1 


0.58 


1.06 


1.14 


Ne nova 3 


V1494 Aql 1999#2 . 


. 1.13 


0.28 


0.1346 


1.24 


0.62 


1.04 


1.12 


Ne nova 3 


V1494 Aql 1999#2 . 


. 1.06 


0.28 


0.1346 


1.22 


0.60 


1.11 


1.26 


Ne nova 2 


V1494 Aql 1999#2 . 


. 0.92 


0.28 


0.1346 


1.17 


0.57 


1.75 


2.00 


CO nova 2 


V2467 Cyg 2007 


. 1.11 


0.35 


0.1596 


1.40 


0.67 


1.08 


1.23 


Ne nova 3 


V2467 Cyg 2007 


. 1.04 


0.35 


0.1596 


1.38 


0.65 


1.16 


1.34 


Ne nova 2 


V2467 Cyg 2007 


. 0.90 


0.35 


0.1596 


1.33 


0.62 


1.67 


1.90 


CO nova 2 


V5116Sgr2005#2 . 


. 1.07 


0.25 


0.1238 


1.15 


0.58 


1.25 


1.43 


Ne nova 3 


V5116Sgr2005#2 . 


. 1.00 


0.25 


0.1238 


1.13 


0.56 


1.30 


1.53 


Ne nova 2 


V5116Sgr2005#2 . 


. 0.85 


0.25 


0.1238 


1.08 


0.53 


2.17 


2.52 


CO nova 2 


V574 Pup 2004 


. 1.05 










1.25 


1.46 


Ne nova 3 


V458 Vul 2007 


. 0.93 


0.74 d 


0.5895 


3.51 


1.40 


1.26 


1.53 


CO nova 4 


V1974 Cyg 1992 


. 1.05 


0.15 


0.0813 


0.85 


0.44 


1.0 


1.17 


Ne nova 2 


V1668 Cyg 1978 


. 0.95 


0.29 


0.1384 


1.21 


0.58 


1.1 


1.3 


CO nova 3 


VI 500 Cyg 1975 


. 1.15 


0.29 


0.1396 


1.28 


0.64 


0.71 


0.8 


Ne nova 2 



a ejected mass by winds 
b envelope mass at the optical maximum 
c see TablefTlfor chemical composition 
upper limit from q = M2/A/WD < <7crit = 0.8 

observational points available before 77 days after the out- 
burst. To understand the characteristic properties of V598 Pup 
light curves, we add observational points of VI 500 Cyg that 
show a very similar light curve in the later phase. We can 
see that two points of V mag in the early phase are just on 
our free-free model light curve but the V data after day ~ 70 
are those in the nebular phase and already deviated from the 
expected free-free light curve. This deviation often occurs 
in novae when V magnitude is contaminated by strong emis- 
sion l ines such as [O III] in the nebular phase (see Read et al. 
200St for emission l ines of V5 98 Pup), as alread y discusse d 
in lHachisu & Kato I d2006l 120071) and lHachisu etaTI (f2008b . 
Thus we regard that the light curve of V598 Pup decays very 



similarly to that of VI 500 Cyg, which follows our universal 
decline law from day ~ 5 to day ~ 70. 

No abundance determinations are available for V598 Pup 
so far. We here assume two typical cases, one is for oxygen- 
neon novae of "Ne nova 2"and the other is for carbon-oxygen 
novae of "CO nova 2" as listed in Table Q] 

The emergence/decay times of supersoft X-ray are a 
good indicator of t he WD mass (lHachisu & Kato I [2006; 
lHachisu et al. 1 12008:). We have calculated several models for 
different WD masses with these two chemical compositions. 
Among the models of Mwd = 1-3, 1.28, and 1.25 Mq, we 
have selected Mwd = 1 -28 M Q as a best-reproducing model 
for "Ne nova 2," but, M W d = 1 -2 M for "CO nova 2" as a 
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best-reproducing one among Mwd = 1-25, 1.2, and 1.15 M©, 
which are not shown in the figure. The best models are se- 
lected by careful eye fitting with the supersoft X-ray data be- 
cause there are only three points. The optical fitting cannot be 
used for the selection because all of three model light curves 
go through the two observational V points as seen in Fig- 
ure [16] We define the "error" in the WD mass determination 
(±0.03 M©) from larger and smaller WD masses, Mwd = 1.3 
and 1 .25 M©, safely cover the decay phase of supersoft X-ray 
flux for the fixed chemical composition. 

The estimated WD mass of 1 .2 M© for "CO nova 2" ex- 
ceeds the upper li mit mass for CO WD s born in a binary, i.e., 
Mco < 1 07 M© (lUmeda et al. II 19991) . so we regard that the 
WD of V598 Pup is very likely an ONeMg WD. Therefore, 
we adopt the 1 .28 M© WD of "Ne nova 2" as a best model, 
which has f bre ak = 40, ? w i n d = 90, and fa-buming = 150 days after 
the outburst (see Figure [Toll. 

We have also calculated the ejecta mass of M w j n( j ~ 3 x 
10~ 6 MfT), which is lost by the optically thick w ind (see, e.g., 
lKato& Hachlsulll994t lHachisu & Katol 120061. for more de- 
tail of ejecta mass calculation). However, it should be noted 
that M w i n( j could be the subject to significant uncertainty be- 
cause we missed the real optical maximum. Various physical 
parameters of novae thus obtained are summarized in Tables 
|6] and Q 

We are not able to determine ?2 and from the light curve 
because of lack of visual magnitudes around the optical peak. 

The mass determination of a WD depends weakly on 
the chemical composition of the envelope, especia l ly, on 
the hydrogen content X (lHachisu & Kato I 120061 12007b 
lHachisu et aini2008t lHachisu & Kato 1 120091) . We have esti- 
mated this dependency from two other chemical composition 
models, i.e., "Ne nova 1" and "Ne nova 3" in TableQ] We sim- 
ilarly obtain a best-reproducing model of Mwd =1-2 Mq for 
X = 0.35 ("Ne nova 1") and M WD = 1 -32 M© for X = 0.65 ("Ne 
nova 3"). The resultant dependency can be approximated as 

Mwd(X) » {1.28 + 0.4(X-0.55)} M , (40) 

for V598 Pup. Here, we have already obtained Mwd(^ = 
0.55) = 1.28 Mq. Therefore, we may conclude that the WD 
mass is Mwd = 1 .28 ±0.04 M© for a typical hydrogen content 
between X = 0.45-0.65 for V598 Pup. The WD mass esti- 
mate depends also on the CNO abundance, but the difference 
among various Xqno models in Table [JJ is smaller than the 
difference among var ious hydrogen content {X) models (see 
Hachi su & Kato 1120061 for more detail of X dependency). 

The distance to V598 Pup can be estimated by two meth- 
ods even if the optical maximum is not available. The first 
one is to use the absolute magnitude of our free-free model 
light curves in Table[3] We use our best-reproducing model in 
Figure [16] and measure the brightness at the end of the free- 
free (FF) light curve, m w , and obtain the distance modulus to 
V598 Pup as 

[(m-M) v ] FF = m w -M w = 11.8-0.1 = 11.7. (41) 

Here we adopt m w =11.8 from Figure [TBI and use the ab- 
solute magnitude of M w = 0.1 at the end point of 1.28 M© 
WD model light curve, which is calculated by interpolation 
between the 1.25 M© and 1.3 M© WD models in Table 
[3] We call this "FF method" as introduced in Section 13.61 
It should be addressed that the FF method is rather robust 
and gives (m-M) v = 11.6-0.0= 11.6 for 1.3 M© WD or 
(m-M) v = 12.2-0.4 = 1 1 .8 for 1 .25 M© WD even if we miss 



a best model. We also emphasize that this FF method does not 
require the optical maximum like the MMRD relations, so it 
can be used even if the optical maximum is not available like 
V598 Pup. 

The next method is to use direct comparison with the other 
nova light curve the distance modulus of which is alre ady 
known. We have already done this method in Section 13.61 
(see Figures [7] [10] and [Till. The direct comparison between 
V1500 Cyg and V598 Pup in Figure [16] gives the difference 
of distance moduli of these two novae, i.e., 

[51og(d/10) + A v ] vl500 Cyg -[51og(rf/10) + A v ] v598 Pup = 0.5. 

(42) 

When we adopt (m-M)v = 12.5 from Equation d20b for 
VI 500 Cyg, we obtain 

[(m-M) v ] hC = [51og(d/10)+A v ] y598 pup = 12.0, (43) 

where LC means the direct light curve (LC) fitting. It should 
be again addressed that the LC method does not need the op- 
tical maximum like the MMRD relations, so it can be used 
even if we missed the optical maximum like V598 Pup. 

The difference of 0.3 mag in the distance modulus esti- 
mation can be understood as follows. These two nova light 
curves look very similar in the later phase but the fitting 
mass is different between 1.28 M© (V598 Pup) and 1.15 M© 
(V1500 Cyg). The absolute magnitude at 20 days, for exam- 
ple, is about 0.4 mag brighter in the 1.15 M© WD model than 
in the 1 .25 M© as shown in Figure [7] We have experienced 
that the same thing happened when we use fitting only of opti- 
cal light curves as shown in Figure[8] In this figure, we obtain 
± 0.4 mag larger (m-M)v when the fitting mass is missed by 
^0.1 M©, respectively. In this sense, the LC method may 
have an error of 0.4 mag or so even if two nova light curves 
look very similar in the later phase. 

The above two estimates are listed in Table[8] being roughly 
consistent with each other, i.e., (m-M)y =11.9 ±0.2 (aver- 
aged with equal weight). These two dista nce moduli give th e 
distance of d = 2.1 ±0.2 kpc for A v = 0.27 dRead et al. 2008). 
This distance is co nsistent with various estimates discussed by 
lReadetal.ld2008l) . 

6.2. V382 Vel 1999 

V382 Velorum (Nova Velorum 199 9) was discovered in- 
depen dently by Williams and Gilmore (Williams & Gilmorel 
1999). The nova reached my =2.7 at maximum on 1999 
May 23 and then decayed to m v = 16.6 (IWoudt et al. 1120051) 
as shown in Figure rrj V382 Vel is a fast nova and identi- 
fied as a neon nova (Woodward et al. 1999). We here adopt 
a set of chemical composition of a typical neon nova, "Ne 
nova 2" in T able [JJ based on the composit i on ana lyses by 
IShore et al. I (I2003F and bvlAugusto & Diaz I (120031) (see Ta- 
ble 1 of Hachis u & Kato 1120061) . We also assume that the out- 
burst day is t 0B =JD 2451319.0 (1999 May 20.5 UT) from 
information of IAU Circulars 7176 and 7184, i.e., m wls > 12 
on May 20.57 UT, m vis = 4.3 on May 21.62 UT, m vis = 3.6 on 
May 21.725 UT, and m v = 3.1 on May 22.396 UT. 

Supersoft X-rays were detected by Bepp oSAX about six 
months after the outburst (lOrio et al. II2002I) . After that, the 
flux decayed rapidly as shown in Fi gure [T7l The Chandra ob- 
servations bv lBurwitz et al. I ([2002) and bv lNess et aTl d2005l) 
suggested that the hydrogen burning ended around 220 days 
after the outburst. The theoretical supersoft X-ray phase is 
very sensitive to the WD mass whereas it depends weakly 
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FIG. 17. — Same as Figure [16] but for V382 Vel. Observational data of 
optical V (boxes with sharp corners) and R (circles with a plus), near infrared 
/ (star marks), and supersoft X-ray (large diamonds). We plot three model 
light curves of Mwd = 1-25 (thick dash-dotted lines), 1.23 (thick solid lines), 
and 1.2 Mq (thin solid lines) WDs for the envelope chemical composition 
of "Ne nova 2." We select the 1.23 Mq WD as a best model among three 
of 1.25, 1.23, and 1.2 Mq WDs and can reproduce the supersoft X-ray data 
(three points). The supersoft X-ray data is sandwiched between the 1.25 and 
1.2 Mq WDs, so we safely regard that M W d = 1.23 ± 0.03 Mq. Observa- 
tional V382 Vel data of V, R, and / magnitudes are taken from IAU Circulars 
7176, 7179, 7196, 7209, 7216, 7226, 7232, 7238, and 7277. The darkest ob - 
servational V magnitudes at day ~ 1850 are taken from Woudt et al. 1 )20051) . 
Visual magnitudes (small circles) are from AAVSO. The X-ray absorbed flux 
dat a (large open d i amond s) of V382 Vel are taken from Olio et al. (2002) 
and Burwitz et al. (2002) without error bars. 



on the chemical composition of the envelope. Our X-ray 
light curve fitting suggests that 1) the WD mass is as massive 
as M WD = 1.23 ±0.03 Mq and 2) the supersoft X-ray phase 
started from day ~ 120. We have also estimated the ejecta 
mass of M wind ~ 4.3 x 10" 6 Mp (Table CD for the 1.23 M Q 
WD model. iDella Valle et al. I (120021) derived the ejecta mass 
of M e j ~6x 1O~ 6 M0 from the flux of Ha emission, which is 
consistent with our estimate. 

Our model of 1.23 Mq WD shows fb rea k = 52, f w ; n d = 120, 
and fH-buming = 220 days after the outburst (see Figure fTTT i. 
Since the V magnitude reached its maximum to = 2.6 days af- 
ter the outburst, we obtain "intrinsic" ti = 6.6 and t$ = 12.4 
days along our model light curve. For this nova, since the ob- 
served data decays almost along our model light curves, the 
"face" values of ?2 and are practically the same as those 
of our "intrinsic" ti and tj,. It should be noted here that the 
assumed epoch of f B = JD 2451319.0 (1999 May 20.5 UT) 
hardly affect our estimation of ?2 and times because we do 
not use the value at ton. Usually we move the model light 
curve "back and forth" and "up and down" against the ob- 
served value as shown in Figure [8] and this process simply 
means that we fit our model light curve with the decay phase 
after maximum. 

In Figure [17] we add V and y light curves of VI 500 Cyg 
for comparison. From the difference between V and y bands, 
we can roughly evaluate the influence of strong emission lines 
that contribute extra flux to V magnitude above free-free con- 
tinuum flux. Comparing the light curves of VI 500 Cyg, we 
regard that V382 Vel light curve is contaminated by strong 
emission lines from ~ 40 days after the outburst; all the V, R, 
and / magnitudes begins to diverge from each other and start 
to deviate from our universal decline law (thick solid line). 

The WD mass estimate depends weakly on the chemical 
composition. Changing hydrogen content X, we obtain a best- 
reproducing model of Mwd =1-13 Mq for X = 0.35 ("Ne nova 



1") and M WD = 1-28 M for X = 0.65 ("Ne nova 3"). Thus, 
we obtain the dependency of the WD mass on X, i.e., 

Mwd(X) w {1.23 + 0.5(^-0.55)} M Q , (44) 

for V382 Vel. This relation on the WD mass is a bit steeper 
than that for V598 Pup. We may also conclude that the WD 
mass is Mwd = 1-23 ±0.05 M© for a typical hydrogen con- 
tent between X = 0.45-0.65 for V382 Vel, which corresponds 
to the div ersity in the composition determinati on between 
X = 66 (IShore et al. 1120031) and X = 0.47 (lAugusto & Diaz I 
120031) . 

The distance modulus to V382 Vel can be calculated in the 
same two ways as those in the previous subsection (V598 
Pup). First, the FF method gives 

[(m-M) y ] FF = m w -M w = 11.9-0.5= 11.4, (45) 

for the 1 .23 M© WD model. Even if we use two other models 
of 1.25 Mq and 1.2 Mq WDs in Figure[l7] this method gives 
a similar result of (m-M) v = 1 1 .7 - 0.4 = 1 1 .3 for 1 .25 Mq, 
ox(m-M) v = 12.2-0.7= 1 1 .5 for 1 .2 M WD. 

Second, the LC method also gives us the absolute bright- 
ness, i.e., 



[51og(d/T0)+A v 



VI 500 Cyg 



[51og(d/T0)+A v 



V382 Vel 



= 1.0. 



With (m-M) v = 12.5 from Equation 
above Equation (|46*] | yields 

[(m-M) v ] LC = [51og(d/10)+A v 



(46) 

for VI 500 Cyg, the 



V382 Vel 



= 11.5. 



(47) 



Kaler- Schmidt's law with tj = 12.4 days and D elia Valle 
& Liv io's law with ti = 6.6 days and my max = 2.6 dGilmore I 
1999) give 



[(m-M) 



V.maxJ MMRD ' 



[51og(£//10)+A V 



RD = 11.6(11.4), 
(48) 

where the value outside parenthesis is calculated from Kaler- 
Schmidt's law of Equation d38l and the value in parenthesis 
is calculated from Delia Valle & Livio's law of Equation d39b . 
These two values are consistent with the FF and LC results. 

These four values give a distance modulus of (m-M)v = 
11.5 ±0.1 for these four different estimations with equal 
weight and a distance of d = 1.5 ±0.1 kpc f or Ay = 3.1 x 
E(B-V) = 3.1 x 0.2 = 0.62 (IShore et al. 112003b . 

These values are consistent with the distance of 1 .66 ±0.11 
kpc derived from the maximum magnitude v s. rate of decline 
(MMR D) relation (IDella Valle et al. 1 120021) . On the other 
hand, IShore et al. I (120031) obtained a bit larger distance of 
d = 2.5 kpc, assuming that the UV flux of V382 Vel is equal to 
that of VI 974 Cyg. They took a distance of 3. 1 kpc to VI 974 
Cyg but this value is much l arger than a reasonable one of 
1.8 kpc (Choc hol et al. Ifl997l) . If we take the distance of 1.8 
kpc instead of 3.1 kpc to V1974 Cyg, Shore et al.'s method 
gives a much shorter distance of 1.5 kpc to V382 Vel, which 
is consistent with our estimates. 

Finally, we introduce some characteristic features concern- 
ing binary nature, that is, the emergence of the companion 
from the WD photosphere. We define this characteristic time 
as the emergence time of the companion, i.e., f em erge- The or- 
bital period of P„ r b = 0.1461 days (3.51 hr) was derived by 
Balm an et al. I ((2006) from the orbital modulations with a full 
amplitude of 0.02-0.03 mag while a bit longer orbita l perio d 
of P 0lh = 0.1581 days was obtained by IWoudt et al. I (120051) . 
So, we estimate the epoch when the companion emerges from 
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FIG. 18. — Same as Figure [16] but for V4743 Sgr. We plot three model 
light curves of Mwd = 1.2 (thick dash-dotted lines), 1.15 (thick solid lines), 
and 1 . 1 Mq (thin solid lines) WDs for the envelope chemical composition 
of "Ne nova 2." We select the 1.15 Mq WD as a best model among three 
of 1.2, 1.15, and 1.1 Mq WDs and can reproduce the supersoft X-ray data. 
Visual (small open circles) and V (filled squares) data of V4743 Sgr are taken 
from IAU Circulars 7975, 7976, and 7982, and from AAVSO. The X-ray flux 
data (large open diamonds) of V4743 Sgr without erro r bars are taken from 
IQrio et al. I f2003T) . lQrio" & Tepedelenlioglu I <2004l) and lNess et aTl I2007al) . 

the nova envelope. If the mass of the donor star can be esti- 
mated from Warner's (1995) empirical formula, i.e., 

Ml « 0.065 ( J°EL V ' , for 1 .3 < < 9 (49) 

Mq \ hours / hours 

we have M2 = 0.31 Mq, which corresponds to the separation 
of a = 1.35 R for M\ = 1.23 Mq, and the effective Roche 
lobe radius of the primary component (WD) of R\ = 0.68 Rq. 
It is about f em erge ~ 32 days when the photospheric radius of 
the nova envelope shrinks to near the orbit (when R p ^ w a). 
This becomes femerge ~ 28 days if we take the orbit plus the 
companion's radius (when R p h w a + R^)- 

Hard X-ray (2-4 keV) flux of V382 Ve l was detected 
by ASCA and RXTE (iMukai & lshida 1 120011) ab out 20 days 
after th e outburst (day 18 with BeppoSAX by lOrio et aLl 
2001b) and reached maximum at day ~ 30 and stayed at 
the same level and th en began to decrease at day ~ 50. 
Mukai & Ishidaj (12001b interpr eted the origin of hard X-rays 
by a model of internal shocks (|Friediung"lll987l) . This shock 
may be for med by collision betwe en two ejecta shells as ex- 
plained by Mukai & Ishida (2001), but there is another pos- 
sibility that the shock is formed by collision between nova 
win ds (optically thick winds) and the companion star (see, 
e.g. JHachisu & Kato ll200i r2006: Hac hisu et al. ll2008h . 

If it is the second case, the hard X-ray emergence should be 
coincident with the emergence of the companion from the WD 
photosphere because hard X-ray is probably absorbed deep 
inside the nova photosphere. The hard X-ray emergence time 
of ~ 20 days is roughly consistent with f em erge ~ 30 days. On 
the other hand, the hard X-ray flux should be declining as 
the wind mass-loss rate is decreasing. This corresponds to 
the break point of free-free emission light curve, because the 
wind mass-loss rate is rapidly decreasing after the break point 
of fbreak = 52 days. Therefore, the decrease in the X-ray flux at 
day 50 is also consistent with our wind model. After fbreak, the 
hard X-ray becomes weak but lasts about 120 days after the 
outburst, in other words, until the optically thick wind stopped 
at fwind = 120 days. 

6.3. V4743 Sgr 2002#3 



V4743 Sgr (Nova Sagittarii 2002 No.3) was discovered by 
Haseda on 2002 September 20 431 U T (JD 2452537.931) 
at mag about 5.0 (lHaseda et all 120021) . iKato et al. I ( 120021) 
reported a prediscovery magnitude of mag 5.5 on Septem- 
ber 18.465 UT. This o bject was not detected by Brown 
dBrown & Pearce 1120021) on September 9.6 UT (limiting mag 
about 12.0), so we here assume that f OB = JD 2452533.0 (2002 
September 15.5 UT) is the outburst day. 

Bright supersoft X-ray ph ases were reported at d ay ~ 1 80 
(2003 March 19 396 UT) by lStarrfield etaTI J2003h (see also 
Ness et al. 1120031) a nd at day ~ 200 (2003 Ap r il 4.92 1 UT) by 
Orio et al. I (120031) . lOrio & Tepedelenlioglu I (120041) reported 
that the supersoft X-ray flux has already declined on 2004 
September 30.783 UT (day - 750) by a factor of 1000. Figure 
rT~8l shows the optical and supersoft X-ray light curves. The 
supersoft X-ray phase started at least ~ 180 days after the 
outburst and ended at day ~ 400. 

The model of 1.15 Mq WD is a best one among three 1.2, 
1.15, and 1.1 Mq WDs, as shown in Figure [181 Here we 
assume a chemical composition of typical neon novae, "Ne 
nova 2." The 1.15 Mq WD model has fb rea k = 70, f w i n d = 180, 
and fH-buming = 380 days. Adopting a magnitude of 5.0 on day 
3.6 along our model light curve of 1.15 Mq WD, we obtain 
?2 = 9.4 days and t$ = 17 days as listed in Tabled We have esti- 
mated the ejecta mass of M w j n d ~7x 10~ 6 Mq . The visual (V) 
magnitudes are roughly fitted with our model light curve until 
about day 40, but are gradually departing from it after that. 
At this stage, the nova probably entered the transition/nebular 
phase and the deviation comes from the contribution of strong 
emission lines. 

We have also checked the accuracy of WD mass determina- 
tion that comes from the ambiguity of chemical composition. 
We obtain a best-reproducing model of Mwd = 1 03 Mq for 
X = 0.35 ("Ne nova 1"), and M WD = 1.21 Mq for X = 0.65 
("Ne nova 3"). Then, we have an approximate relation of 

M WD (X) w {1.15 + 0.6(X-0.55)} M , (50) 

for V4743 Sgr. This relation of the WD mass is steeper (on 
X) than that for V598 Pup and V382 Vel. We may also con- 
clude that the WD mass is M WD = 1 . 15 ±O.O6M for a typical 
hydrogen content between X = 0.45 - 0.65 of V4743 Sgr. 

From the FF method, we obtain the distance modulus to 
V4743 Sgr, i.e., 

[(m-M) v ] FF = m vl -M vl = 14.7-0.9 = 13.8, (51) 

for the 1.15 Mq WD. Two other attendant models of 1 .2 Mq 
and 1 . 1 Mq WDs in Figure [18] give similar results of (m - 
M) v = 14.3-0.7= 13.6 for 1.2 Mq or (m-M) v = 15.2-1.2 = 
14.0 for 1.1 Mq WD. 
The LC method gives 

[51og(d/10) + A v ] vl500Cyg -[51og(d/10)+A v ] V4743 Sgr = -l.l 

(52) 

With (m-M) v = 12.5 from Equation (HI for V1500 Cyg, we 
obtain 

[(m-M) v ] LC = [51og(d/10)+A v ] v4743 Sgr = 13.6. (53) 

Kaler-Schmidt's law with f3 = 17 days and Delia Valle & 
Livio's law with f2 = 9.4 days and my max = 5.0 yield 

[(^-M) V!raax ] MMRD =[51og(,//10) + A v ] MMRD =13.7(13.7), 

(54) 
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FIG. 19. — Same as Figure [T6l but for V1281 Sco. Observational data of 
optical V (red filled squares) are taken from AAVSO and Association Fran- 
caise des Observateurs d'Etoiles Variables (AFOEV). The X-ray count rate 
data (large open diamonds) with error bars (magenta bar s) of VI 281 Sco ar e 
taken from an automatic analyzer of the Swift web page I Evans et al. 2009). 
We plot three model light curves of Mwd = 115 (thick dash-dotted lines), 
1.13 (thick solid lines), and 1.1 Mq (thin solid lines) WDs for the envelope 
chemical composition of "Ne nova 2." We select the 1.13 Mq WD as a best 
model among three of 1.15, 1.13, and 1.1 Mq WDs and can reproduce the 
supersoft X-ray data. 

being consistent with the FF and LC results. Then, the dis- 
tance modulus to V4743 Sgr is (m-M) v = 13.7 ±0.1 and 
the distance is d = 3.8 ±0.2 kpc for A v = 2>.\E(B-V) = 
3.1 x 0.25 = 0.78 dVanlandingham et al. Il2007l) . Our estimate 
is also consistent with th e recent results of E(B - V) = 0.25 
and d = 3.9 ± 0.3 kpc bv lVanlandingham et al. I (120071) . 

The o rbital period of P nr h = 0.2799 days (6.72 hr) was de- 
rived by Kang et al. ( 2006) from the orbital modulations with 
a full amplitude of 0.05 mag (2003) - 0. 15 mag (2005). Then, 
we estimate the epoch when the companion emerges from the 
nova envelope. From Equation j49l , we have M2 = 0.70 Mq, 
which corresponds to a = 2.2 Rq and R^ = 0.93 Rq for Mi = 
1.15 Mr. It is r eme 



36 days at R p h ~ a. 



6.4. V 1281 Sco 2007 #2 



VI 281 Sco (Nova Scorpii 2007 No. 2) was discovered by 
Nakamura on 200 7 February 19.86 UT (JD 2454151.3593) 
at mag about 9.3 ( Yama oka et al. 1120071). This object was 
not detected by Nakamura (lYamaoka et al. 1120071) on Febru- 
ary 14.86 UT (limitin g mag about 12.0) and also by Fu- 
jita (iNaito et al. 112007b on February 18.85 UT (limiting mag 
about 1 1 .6), so we here assume that f B = JD 2454148.5 (2007 
February 17.0 UT) is the outburst day. 

A bright supersoft X-ray phase was repor t ed at d ay ~ 
340 (2008 January 24.18 UT) by iNess et al. I d2008ah . To 
search for the end epoch of supersoft X-ray phase, i.e., 
^H-huming, we use an automatic analyzer in the Swift web page 1 
(Evan s et al. 120091) . The X-ray (0.3-10keV) countrates with 
the Swift XRT are plotted in Figure [19] together with its op- 
tical light curve. We have checked the early data points of 
X-ray obtained b y automatic analyzer are very similar to the 
values reported bv lNess et al. I(l2008al) . The epoch of ?H-buming 
is about 425 days after the outburst, if we regard that the last 
observation indicates the decay of supersoft X-ray flux. 

Our best-reproducing light curve of a 1.13 M© WD is plot- 
ted in Figure[l9]for a chemical composition of a typical neon 
nova, "Ne nova 2", which shows fb rea k = 78, f w ; n d = 190, and 

1 http://www.swift.ac.uk/ 



fH-buming = 425 days after the outburst. We determined the WD 
mass to be 1.13 ±0.03 M© mainly from the supersoft X-ray 
data. Adopting a magnitude of 9.2 on day 3.0 along our model 
light curve of 1.13 M© WD, we obtain f 2 = 9.5 and £3 = 16.7 
days as listed in Table|6] We have estimated the ejecta mass of 
^wind ~9x 10~ 6 M©. The visual (V) magnitudes are roughly 
fitted with our model light curve until about day 80, but are 
gradually departing from it after that. This deviation comes 
from the contribution of strong emission lines. 

We have also examined the accuracy of WD mass determi- 
nation that comes from the ambiguity of chemical composi- 
tion. We obtain a best-reproducing model of Mwd = 1-01 M© 
for X = 0.35 ("Ne nova 1"), and M WD =1.2 M© for X = 0.65 
("Ne nova 3"). Then, we have 



M WD {X) « {1. 13 + 0.6(^-0.55)} M G 



(55) 



for VI 281 Sco. We may conclude that the WD mass is Mwd = 
1.13 ±0.06 Mq for a typical hydrogen content between X = 
0.45 -0.65 of V1281 Sco. 

From the FF method, we obtain the distance modulus to 
V1281 Sco, i.e., 



[(m-M) v ] FP = m vl -M w = 18.8-1.0= 17.8, 



(56) 



for the 1.13 M© WD model. Two other attendant models of 
1.15 M© and 1.1 M© WDs in Figure [19] give similar results 
of (m-M) v = 18.6-0.9 = 17.7 for 1.15 M© or (m-M) v = 
19.1-1.2= 17.9for 1.1 M© WD. 
The LC method gives 



[5\og(d/10)+A v ] V150Q Cyg - [51og(<//10)+A v 



I V1281 Sco 



-5.0. 
(57) 

With (m-M) v = 12.5 from Equation (E3 for V1500 Cyg, we 
obtain 



[(m-M) v ] LC = [51og( fl f/10)+A v 



V1281 Sco 



= 17.5. 



(58) 



Kaler- Schmidt's law with £3 = 16.7 days and Delia Valle & 
Livio's law with t^ = 9.5 days and my. max = 9.2 (at frj = 3.0 
days) give 

[(m-M) v , max ] MMRD = [51og(rf/10)+A v ] MMRD = 17.9 (17.9), 

(59) 

being consistent with the FF and LC results, where the differ- 
ence between the LC method and the others is about ± 0.4 
ma g but within the ambiguity in this method (see Sections l4~2l 
and !6. It . Then the distance modulus is (m-M)y = 17.8 ±0.2 
and the distance to VI 281 Sc o is d = 13 ± 1 kpc for Ay = 
3AE(B-V) = 3.l x 0.7 = 2.17 dRussell et al. Il2007al) . 

6.5. V597 Pup 2007#1 

V597 Pup (Nova Puppis 2007 No.l) was discovered by 
Pereira on 2007 November 14, 23 UT (JD 2454418.73) at mag 
about 7.0 (Perei raet al. 11200 7). This object was not detected 
dPereiraet all 120071) on November 6.23, 7.22, 8.23, 10.23, 
and 11.22 UT (limiting mag about 8), so we assume that 
? 0B = JD 2454416.0 (2007 November 1 1.5 UT) is the outburst 
day. 

A bri ght supersoft X-ray phase was reported by Nes s et al. I 
( 2008c). To search for the end epoch of supersoft X-ray phase, 
i.e., f H-hiirning, we use an automatic analyzer in the Swift web 
page tevans etal. 1 12009). The X-ray (0.3 - 10 keV) count 
rates with the Swift XRT are plotted in Figures l20l and |2TI as 
well as the optical light curves. The X-ray data points show 
that the turnoff of supersoft X-ray is fH-buming ~ 430 - 460 
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FIG. 20. — Same as Figure[l6] but for V597 Pup. We plot three model 
light curves of Mwd = 1.1 (thick dash-dotted lines), 1.05 (thick solid lines), 
and 1.0 Mq (thin solid lines) WDs for the envelope chemical composition 
of "CO nova 2." We do not select the best model among these three of 1.1, 
1.05, and 1.0 Mq WDs, because the supersoft X-ray duration is too short 
to be comparable with the observation. The X-ray count rate data (0.3 — 10 
keV: large open diamonds) with error bars (mage nta bars) ar e taken from an 
automatic analyzer of the Swift web page ( Evans et al. 2009). Observational 
V597 Pup data of V (filled squares) magnitudes are taken from AFOEV and 
AAVSO. 
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FIG. 21. — Same as Figure [20l but we plot three model light curves of 
Mwd = 1.25 (thick dash-dotted lines), 1.2 (thick solid lines), and 1.15 Mq 
(thin solid lines) WDs for the envelope chemical composition of "Ne nova 
3." We select the 1.2 Mq WD as a best model among three of 1.25, 1.2, and 
1.15 Mq WDs and can reproduce the supersoft X-ray data. 

days and the turn-on is probably f w j n d ~ 150 - 170 days. This 
is a rather rare case in which both the turn-on and turnoff 
times of supersoft X-ray are specified. 

The chemical composition of ejecta is not known, so we 
first assumed a chemical composition of "CO nova 2" and 
made model light curves as shown in Figure|20] The durations 
of supersoft X-ray phases are too short to be comparable with 
the observation. This is because the duration of supersoft X- 
ray is shorter for a smaller hydrogen content of X. 

Our reasonably-reproducing model light curve of a 1 .2 Mq 
WD is plotted in Figure|2T|for a chemical composition of "Ne 
nova 3", which shows f brea k = 74, f wind = 200, and fH-buming = 
430 days after the outburst. We estimated the WD mass to be 
1 .2 ± 0.05 Mq mainly from the supersoft X-ray data. Adopt- 
ing a magnitude of 8.4 on day 2.7 as the maximum magnitude 
along our model light curve of 1 .2 Mq WD, we obtain ti = 8.7 
and ?3 = 16.5 days as listed in Table We have estimated the 



ejecta mass of M w ; n d ~ 1 x 1O~ 5 M . The optical light curve 
shape around the peak deviates largely from our model light 
curves and is similar to that of V1500 Cyg. It is highly likely 
that V597 Pup is a superbright nova. After the superbright 
phase, V magnitudes are roughly fitted with our model light 
curve until day ~ 100. 

We have also checked the accuracy of WD mass determina- 
tion from the ambiguity of chemical composition. For chem- 
ical compositions of "Ne nova 1" and "Ne nova 2," however, 
we could not obtain reasonable fits like those as shown in Fig- 
ure |2T|because supersoft X-ray durations of these models are 
too short to fit with the supersoft X-ray observation like in 
Figure [20] Therefore, we may conclude that the WD mass is 
Mwd = 1.2 ±0.05 Mq and that the composition is relatively 
hydrogen-rich like X ~ 0.65. 

From the FF method, we obtain the distance modulus to 
V597 Pup, i.e., 

[(m-MV] FF = m w -M w = 17.8-0.9= 16.9, (60) 

for the 1.15 Mq WD ("Ne nova 2") model. Two other at- 
tendant models of 1.2 Mq and 1.1 Mq WDs (not shown 
in the figure for "Ne nova 2" models), give similar results 
of (m-M)v = 17.5-0.7 = 16.8 for 1.2 Mq or (m-M) v = 
18.3-1.2 = 17.1 for 1.1 Mq WD. From Figure GQ] we also 
have (m-A0v = 17.3-0.5 = 16.8 for 1.1 Mq, (tn-M) v = 
17.8-0.8 = 17.0 for 1.05 M , or (m-M) v = 18.3-1.1 = 17.2 
for 1 .0 Mq WD with "CO nova 2." 
The LC method gives 



[51og(£//10)+A v 



I VI 500 Cyg 



[51og(<i/10) + A v ] V597 pup = -4.2. 

(61) 

With (m-M) v = 12.5 from Equation (|20j for V1500 Cyg, we 
obtain 



[(m-M) v ] LC = [51og( t //10)+A v ] V597 pup = 16.7. 



(62) 



Kaler- Schmidt's law with £3 = 16.5 days and Delia Valle & 
Livio's law with t<i = 8.7 days, and both with my. max = 8.4 (at 
to = 2.7 days) yield 



[{m-M) Vr 



Jmmrd 



= [51og(<//10)+A v ] 



MMRD 



= 17.1(17.1), 
(63) 

being consistent with the FF and LC results. Then we have 
(m-M) v = 17.0 ± 0.2 and the distance to V597 P up,d = 16 ± 
2 kpc for A v = 3AE(B-V) = 3.1 x 0.3 = 0.93 dNess et al. I 
20083). 

The orbi tal period of P nr b = 0.11 12 days (2.6687 hr) was 
derived by Warner & Woudt (2009) from the orbital modu- 
lations with a full amplitude of 0.2 mag (2008) - 0.6 mag 
(2009). Then, we estimate the epoch when the companion 
emerges from the nova envelope. From Equation (|49l , we 
have M2 = 0.22 Mq, which corresponds to a = 1.1 Rq and 
R* = 0.58 ^ for Mi = 1.2 Mq. It is f e mer g e ~ 52 days at 

War ner & Woudt I (120091) reported that a repetitive hump 
with an amplitude of ~ 0.2 mag is seen but no obvious eclipse 
features in 2008 (~ 125 days after outburst). Their Fourier 
analysis of light curves provide an orbital period of 2.67 hr. 
After that, their 2009 observation revealed broad eclipses with 
~ 0.6 mag deep (~ 470-500 days after outburst). They in- 
terpreted that the drop in brightness (2.5 mag) between these 
two observations in 2008 and 2009 was due largely or entirely 
to the dispersal and cooling of nova ejecta far from the cen- 
tral binary and the eclipse and their related humps could well 
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FIG. 22. — Optical and supersoft X-ray light curves of V1494 Aql. Opti- 
cal bands of visual (small open circles), V (filled squares), R (circles with a 
plus), and near IR band of / (small open diamonds) magnitudes are shown, 
which are taken from AAVSO and VSNET. The X- ray count rate data (larg e 
open diamonds) without error bars are taken from Rohrbach et al. 1 120091) . 
We plot three model light curves of Mwd = 108 (thick dash-dotted lines), 
1.06 (thick solid lines), and 1.04 Mq (thin solid lines) WDs for the envelope 
chemical composition of "Ne nova 2." We select the 1.06 Mq WD as a best 
model among three of 1.08, 1.06, and 1.04 Mq WDs and can reproduce the 
supersoft X-ray data. Observational V magnitud es of V197 4 Cyg are added 
for comparison, which are taken from Chochol et al. I j!993l) . 

have been present even in the earlier 2008 light curves, that 
is, there are some features resembling such diluted structures. 
This interpretation is consistent with the emergence of com- 
panion on ~ 50 (Rph = a) or ~ 60 days (R p h = a-R.2) after 
outburst. Moreover, the WD photosphere is larger than that of 
the companion, i.e., /? p h > R2, during the optically thick wind 
phase, at least, until 130 days after outburst. After the winds 
stop on day ~ 200, the WD photosphere rapidly shrinks to- 
ward its origina l size, that is, it becomes a half by one month 
(see, e.g. JHachisu & Kato 1120031) . 

6.6. V '1494 Aql 1999#2 

V1494 Aql (Nova Aquilae 1999 No.2) was discovered by 
Pereira on 1999 December 1 785 U T (JD 2451514.285) at 
mag about 6.0 dPereira et al. 1 1 19991) . This object was not 
detected by Liller (Fuiii et al. 111999b on November 25.035 
UT (limiting mag about 10.5), so we assume that foB = JD 
245 15 13.0 (1999 November 30.5 UT) is the outburst day. The 
star reached its maximum of m v ; s ss 3.9 on 1999 December 3.4 
UT (JD 2451515.9) about 3 days after the outburst. Then it 
declined by 2 mag in 6. 6 days and by 3 mag in 16 days (e.g., 
iKiss & Thomson 1120001) . The transition oscillation began on 
day ~ 20 and ended on day ~ 130, during which the star ex- 
hibited quasi-periodic oscillations with periods of 7 to 22 days 
and with an amplitude reaching 1.5 mag. 

We plot the optical and supersoft X-ray light curves in 
Figure [22] A bright supersoft X-ray phase was detected 
with Chandra on day 250 (2000 August 6 UT), day 302 
2000 September 2 8 UT), and day 305 ( 2000 October 1 UT) 
Drake et al.ll2003l) . iNess et al. I d2007al) reported that the su- 
persoft X-ray flux had already declined by a factor of 40 on 
day - 730 (2001 November 28 UT). Then the supersoft X- 
ray phase started at least on day ~ 250 and ended before day 

~ 700. 

Iiiim a & Esenoglul d2003l) analyzed their spectra of V 1494 
Aql and estimated the helium/hydrogen ratio to be He/H= 
0.13 ±0.01 by number ratio, which corresponds to Y/X = 
0.52 ±0.04. This helium/hydrogen ratio is close to our chem- 



ical composition of "N e nova 2," in which Y/X w 0.54. 
lArkhipova etaLl d2002l) derived He/H= 0.126, N/H= 0.04 (or 
0.0015), 0/H= 0.002, and Fe/H= 5 x 10" 5 by number ratio. 
This corresponds to X = 0.48, Y = 0.24, X CN o = 0.26, and 
Z = 0.02, (or to X = 0.62, Y = 0.34, X CN o = 0.02, and Z = 0.02). 
We do not know which is better, so that we adopt an arithmetic 
average of these two, i.e., X = 0.55, Y = 0.29, X CN o = 0. 14, and 
Z = 0.02, which is close to "Ne nova 2." Therefore, we assume 
"Ne nova 2"as the envelope chemical composition of V1494 
Aql. 

Our best-reproducing light curve is a 1 .06 Mq WD model 
among 1 .04, 1 .06, and 1 .08 Mq WDs mainly from the super- 
soft X-ray data (Figure[22]), i.e., M WD = 1 .06 ± 0.02 Mq. For 
the 1.06 Mq WD model, we derive fb rea k = 102, f w ; n( j = 260, 
and f H-burning = 670 days. Adopting a magnitude of mv, max = 
5.0 on day 2.7 as the optical maximum along our model light 
curve of 1 .06 Mq WD, we obtain ti = 1 1 .7 days and = 22.1 
days (Table[6]). The visual (V) magnitudes are not well fitted 
with our model light curve, from day 10 to day 150, mainly 
because the light curve has wavy structures and seems to be 
strongly contaminated by emission lines. Our model light 
curve almost follows the bottom line of V magnitudes during 
the transition oscillations, whic h is supported by Figu re [TOl 
for GK Per (see also Figure 2 of Hachisu & Kato 2007). The 
nova entered a nebular phase ^100 days after the outburst and 
the visual magnitudes deviated much from our model light 
curve due mainly to stron g emission lines such as [O III] (e.g., 
Iiiima & Esenoglu 2003). 

We have also estimated the ejecta mass to be M w j n d ~ 1.1 x 
1O~ 5 M (Table [7J, which is lost by the optically thick w ind 
(see, e.g., Kato & Hachis ulflS^ [fi achisu & Katoll2006L for 
more detail), [iiiima & Esenoglul d2003l) derived the ejecta 
mass of M ej ~ (6.2 ± 1.4) x 10" 6 M Q (J/1.6kpc) 2 , being con- 
sistent with our model estimation if we take the distance of 
d = 2.2 kpc as derived below. 

Based on submillimeter- and centimeter-band fluxes, a sin- 
gle M ERLIN image, and optical spectroscopy, Evres et al. 
(2005) concluded that optical spectroscopy indicates contin- 
ued mass ejection for over 195 days, with the material becom- 
ing optically thin sometime between 195 and 285 days after 
the outburst. This is consistent with f w ; n( i = 260 days of our 
1 . 06 Mq WD model. 

iMazuk et al71 ( 120001) reported that VI 494 Aql had entered 
into a strong coronal phase at least 230 days after the outburst. 
This is closely related to the epoch when the optically thick 
winds stop, because the photosphere rapidly shrinks and the 
photospheric temperature quickly increases to ~ 3 x 10 5 K or 
more. 

We have checked the dependency of the WD mass de- 
terminatio n on the hydr ogen content of the WD enve- 
lope. iKamath et al. I d2005l) reported a different result of he- 
lium/hydrogen ratio of the ejecta from He II 4686/H/? ratio, 
that is, He/H= 0.24 ± 0.06 by number ratio, which corre- 
sponds to Y/X = 0.96 ±0.24 and is close to our "CO nova 
2" (see Table Q]). When we adopt a chemical composition of 
"CO nova 2," we obtain a best-reproducing model of 0.92 Mq 
WD. The other physical parameters are also listed in Tables|6] 
and|7] On the other hand, if we assume a chemical composi- 
tion of "Ne nova 3," we obtain a best model of 1.13 Mq WD. 
Finally we have 

Mwd(X) w {1.06 + 0.7(X-0.55)} M , (64) 

for V1494 Aql. Here, we use M WD (Z = 0.65) =1.13 Mq, 
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M WD (X = 0.55) = 1.06 Mq, and M WD (X = 0.35) = 0.92 M . 
Then, we may conclude that the WD mass is Mwd = 1 06 ± 
0.07 Mq for a typical hydrogen content between X = 0.45 - 
0.65. If the hydrogen content is X ~ 0.55, then the WD has 
a mass of 1.06 ±0.03 M and is probably a ONeMg WD. 
On the other hand, when the hydrogen content is as low as 
X ~ 0.35, then the WD is less massive (0.92 ± 0.03 Mq) and 
probably a CO WD. 

From the FF method, we obtain the distance modulus to 
V1494 Aql,i.e., 

[(m-M) v ] FF = m w -M w = 14.9-1.5= 13.4, (65) 

for the 1.06 Mq WD ("Ne nova 2") model. Two other at- 
tendant models of 1 .08 Mq and 1 .04 Mq WDs give simi- 
lar results of (m-M) v = 14.6-1.1 = 13.5 for 1.08 M or 
(m-M) v = 15.0-1.7= 13.3 for 1 .04 Mq WD. 
The LC method gives 

[51og(d/10)+A y ] vl974 Cyg -[51og(rf/10)+Av] vl494Aql = -l.: 

(66) 

With (m-V) v = 12.3 from Equation (O for V1974 Cyg, we 
obtain 

[(m-M) v ] hC = [5log{d/l0)+A v ] vim Aq] = 13.6. (67) 

Kaler- Schmidt's law with ?3 = 22.1 days and Delia Valle & 
Livio's law with t2 = H-7 days and my max = 5.0 (at ?o = 2.7 
days) yield 

[(m-M) v . max ] MMRD = [51og(<//10)+A v ] MMRD = 13.4(13.6), 

(68) 

being consistent with the FF and LC results. Then, we obtain 
the distance modulus of (m-M) v = 13.5 ±0. 1 and the distance 
of d = 2.2 ±0.2 kpc to V1494 Aql for A v = 3.1E(B-V) = 
3.1 x 0.6 = 1.83 (llijima & Esenogl u 2003). This value is 
a bit larger than the dista nce of 1.6 ±0.2 kpc estimated by 
llijima & Esenoglu I (120031) from the MMRD relation. 

The orbi tal period of P nr h = 0.1346 da ys (3.23 hr) was 
derived by iBarsukova & Goranskiil (120031) from the orbital 
modulations with an eclipse depth of ~ 0.5 mag (see also 
iRetter et al.1l2000t iBos et al. 1 120011) . Then, we estimate the 
epoch when the companion emerges from the nova envelope. 
From Equation j49l , we have M2 = 0.28 M@, which corre- 
sponds to a = 1.22 R Q and R\ = 0.60 R Q for M x = 1.06 M Q . 
These give f e merge ~ 72 days at R p h ~ a. Various physical pa- 
ra meters are summari z ed in T ables |6l and |?J 

llijima & Esenoglu I (120031) and [Kamath e faXI (120051) re- 
ported that higher ionization lines became prominent about 
65 days after the optical maximum (i.e., day ~ 70 after the 
outburst). At that time, the photospheric temperature of our 
1 .06 Mq WD model has already increased to T p h ~ 60000 - 
70000 K and the higher ionization lines mentioned above are 
possibly related to this high photospheric temperatures. An- 
other possibility is high temperature optically thin plasma that 
emits hard X-ray mentioned below, because high tempera- 
ture plasma that emits hard X-ray emerges around t emel <, e ~ 70 
days. 

Hard X-ray was dete cted with Chandra long before the su- 
persoft X-ray appeared. Rohr bach et al. I (120091) reported their 
six Chandra X-ray spectra. The first three observations were 
taken with ACIS-I on 134, 187, and 248 days after outburst. 
The count rates are 1.0, 0.69 and 0.53 counts s -1 , respectively. 
(We do not plot these point in Figure |22]because the detector 
ACIS-I is different from that of supersoft X-ray observation, 
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FIG. 23. — Same as Figure l22l but for V2467 Cyg, i.e., optical and super- 
soft X-ray light curves of V2467 Cyg. Optical V (open red squares), visual 
(small red filled squares), R (ocher open circles with plus), and I (magenta 
open diamonds) are taken from IAU Circular 8821, Tomov et al. (2007), and 
AAVSO, while the supersoft X-ray dat a (large black diamo nds) with error 
bars (magenta bars) are compiled from Ness et al. (2008b) and also taken 
from an automatic analyzer of the Swift web page (Evans et al. 2009). We 
plot three model light curves of Mwd = 108 (thick dash-dotted lines), 1 .04 
(thick solid lines), and 1.0 Mq (thin solid lines) WDs for the envelope chem- 
ical composition of "Ne nova 2." We select the 1 .04 Mq WD as a best model 
among three of 1 .08, 1 .04, and 1 .0 M© WDs and can reproduce the supersoft 
X-ray data. 

HRC-S.) They found no significant variations in any of these 
observations. The spectra are dominated by emission lines 
originating from N, O, and Ne. They fit isothermal APEC 
models with the spectra and found best fits with kT = 0.65 keV 
and N u = 3.3 x 10 21 c m 2 (E(B-V) = N H /4.8 x 10 21 = 0.68 
by iNess et aLll2007ah . They reported that, in all cases the 
elemental abundances of O and N had to be significantly in- 
creased. They further reported that on day 248 a bright soft 
component appeared in addition to the fading emission lines. 
When the harder X-ray declined, that is, when the optically 
thick wind stopped, the soft component of X-ray emerged. 
This is very consistent with our collision model between the 
wind and the companion as an origin of high temperature 
plasma that emits hard X-ray. In our model, the hard compo- 
nent of X-rays emerges around f em erge ~ 70 days and attains its 
peak at fb re ak ^100 days and then declines toward f w ; n d ~ 260 
days. 

6.7. V2467 Cyg 2007 

V2467 Cyg (Nova Cygni 2007) was discovered by Tago 
on 2007 March 15.78 7 UT (JD 2454175.287) at 7.4 mag 
dNakano et al. Il20~07bl) . Tago also reported that the nova was 
fainter than 12th magnitude three days ago (JD 2454172.296). 
Here we assume the outburst day of foB = JD 2454172.0 
(2007 March 12.5 UT). In Figure |23] we plot optical and 
near infrared magn itudes taken from IAU Circular 8821, 
iTomov et al. I d2007l) , and AAVSO. The nova reached its max- 
imum of my, m ax « 7.4 at t w JD 2454175.2 (March 16.7 UT). 
Except for some early, low amplitude oscillations, the light 
curve declined monotonically until day ~ 30. Then V2467 
Cyg started its transition phase until day ~ 150 in the form of 
a series of several oscillations with a quasi-period of 20-22 
d ays and a 0.5-1.1 m ag amplitude. 

INess et aL~l (|2008b) obtained two Swift observations. The 
first observation was taken on 2008 June 15.17, about 460 
days after the outburst, and obtained an X-ray count rate of 
0.017 counts s" 1 . About 100 days later (2008 September 23.3 
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UT = 566 days after the outburst), their second observation 
yields a count rate of 0.048 counts s" 1 . The both X-ray spec- 
tra are soft with a blackbody temperature of kT = 34 eV and 
kT = 46 eV, respectively. These two count rates are plotted 
in Figure [23] To search for the end epoch of supersoft X- 
ray phase, i.e., frfcjanninai we use an automatic analyzer in the 
Swift web page dEvans et al. Il2009h . The three X-ray (0.3-10 
keV) data points with the Swift XRT are plotted in Figure l23l 
early tw o of which are the same data as those by Ness et al. 
(2008b). We have checked that the early two data points of 
the automa ti c analy zer are very similar to those obtained by 
iNess et al~l d2008bi) . The epoch of fH-buming is around 760 
days after the outburst. 

Our free-free model light curves are plotted in Figurel23lfor 
a chemical composition of "Ne nova 2." The bottom line of R 
and / magnitudes are well fitted with our model light curve un- 
til about day ~ 100. On the other hand, V magnitude is finally 
departing from the R and / magnitudes from day ~ 150 due 
to strong emission line contributions. The best-reproducing 
WD mass is 1.04 M among 1.08, 1.04, and 1.0 M Q . For the 
1 .04 M Q WD model, we obtain f break w 1 13, ? wind ~ 290, and 
?H-buming ~ 765 days. We also obtain "intrinsic" t2 ~ 12 and 
?3 w 23 days along our model light curve. Here we assume 
that the model light curve reached its maximum (7.9 mag) 
at to ~ 3.3 days after the outburst (at the first observational 
p oint). 

ISwierczvnski et al~l (120081) reported the possible orbital pe- 
riod of V2467 Cyg, P orb = 0.1596 days (3.82 hr), from their 
unfiltered, R, and / band CCD photometry. Similar periods of 
3.57 hr and 3.85 hr had be en already observed in 2007, several 
weeks after maximum, by lTomov et al. I (12007b . We estimate 
the epoch when the companion emerges from the nova enve- 
lope. With the donor mass of M2 = 0.35 Mq estimated from 
Equation ([49), we obtain a = 1.38 R Q and R* = 0.65 R Q for 
Mi = 1 .04 Mq, giving ? e mei-ge ~ 75 days at R p h ~ a- 

We have checked the accuracy of WD mass determination 
from the ambiguity of chemical composition. We obtain a 
best-reproducing model of M W d = 0.90 Mq forZ = 0.35 ("CO 
nova 2"), and M WD = 1.11 Mq for X = 0.65 ("Ne nova 3"). 
Then, we have 

Mwd(X) » {1.04 + 0.7(X-0.55)} M , (69) 

for V2467 Cyg. Thus, we may conclude that the WD mass 
is Mwd = 1 04 ± 0.07 Mq for a typical hydrogen content be- 
tween X = 0.45 - 0.65. If the hydrogen content is X ~ 0.55, 
then the WD has a mass of 1.04 ±0.03 Mq and is proba- 
bly a ONeMg WD. On the other hand, when the hydrogen 
content is as low as X ~ 0.35, then the WD is less massive 
(0.90 ± 0.03 M ) and probably a CO WD. 

From the FF method, we obtain the distance modulus to 
V2467 Cyg, i.e., 

[(m-M) v ] pp = m w -M vl = 18.0-1.7= 16.3, (70) 

for the 1 .04 M Q WD ("Ne nova 2") model. Two other at- 
tendant models of 1.08 Mq and 1.0 Mq WDs give simi- 
lar results of (m-M) v = 17.6- 1.4 = 16.2 for 1.08 M or 
(m-M) v = 18.5-2.0= 16.5 for 1.0 M© WD. 
The LC method gives 

[51og(rf/10)+A v ] vl974 Cyg - [51og(d/10)+Av] V2467 Cyg = -4. 

(71) 

With (m-V)v = 12.3 from Equation dT6} for V1974 Cyg, we 
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FIG. 24.— Same as Figure [22] but for V5116 Sgr. Optical V (filled red 
squares), visual (small red open circles), and supersoft X-ray (large dia- 
monds) light curves of V51 16 Sgr . The discovery magnitude is denoted by 
a large open square (Liller 2005). The V and visual data are taken from 
AAVSO. The X-ray flux data with error bars (ma genta bars) ar e taken from 
an automatic analyzer of the Swift web page l Evans et al. 2009). We plot 
three model light curves of Mwd = 11 (thick dash-dotted lines), 1.07 (thick 
solid lines), and 1.05 Mq (thin solid lines) WDs for the envelope chemical 
composition of "Ne nova 3." We select the 1.07 Mq WD as a best model 
among three of 1.1, 1.07, and 1.05 Mq WDs and can reproduce the supersoft 
X-ray data. 

have 

[(m-M) v ] hC = [51og(rf/10)+A v ] V2467 Cyg = 16.4. (72) 

Kaler-Schmidt's law with £3 = 23.0 days and Delia Valle & 
Livio's law with f2 = 12.0 days and my. max = 7.9 at maximum 
(at to = 3.3 days) give 

[(m-M) v , max ] MMRD = [51og(<i/10) + Av,] MMRD = 16.2 (16.5), 

(73) 

being consistent with the FF and LC results. Then, we obtain 
the distance modulus of (m-M)v = 16.4±0.2 a nd the distance 
ofd = 2.2 ±0.2 kpc to V2467 Cyg for Ay =4.65 dMazuk et al. I 
l20Q7h . 

6.8. V5116Sgr2005#2 

V5116 Sgr (Nova Sagittarii 2005 No. 2) was discovered 
by L iller on July 4 .049 UT (JD 2453555.549) at mag about 
8.0 dLiller 1 120051) and became a supersoft X-ray source 
dSala et al. 1120081) We have al ready published analysi s of op- 
tical light curves for this nova (Hachisu & Kato 2007) before 
the supersoft X-ray was detected by S ala et al. I (|2007). This 
o bject was not detec t ed on June 12 (limiting mag about 1 1.0), 
so lHachisu & Kato I d2007l) assumed that f 0B =JD 2453552.0 
(July 1.5 UT) is the outburst day. Their best-reproducing 
model is a 0.9 M Q WD for the chemical composition of "CO 
nova 2." Adopting 8.0 mag on JD 2453555.549 observed by 
Liller (Liller 2005) as the maximum visual magnitude for 
V51 16 Sgr, Hachisu & Kato derived ?2 = 20 days and £3 = 33 
days. From their best light curve of 0.9 Mq WD, Hachisu 
& Kato predicted that the supersoft X-ray phase starts on 
fwind ~ 310 days and ends on fH-buming ~ 770 days after the 
outburst. 

We plot the supersoft X-ray light curves as well as opti- 
cal in Figure [24] A bright s upersoft X-ray ph ase was first 
detected with XMM Newton dSala et al. 1 12007b on day 674 
(2007 March 5 UT). ISala et al. I d2008l) reported that the X- 
ray light curve shows abrupt decreases and increases of the 
flux by a factor of ~ 8, which is consistent with a periodic- 
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ity of 2 .97 hr, the orbital period suggested by Dobrotka et al. 
(2008), and speculated that the X-ray light curve may result 
from a partial coverage b y an asymm e tric acc retion disk in 
a high-inclination system. iNess et al71 (l2007d) also reported 
the supersoft X-ray de tection with Swift on day 768 (2007 
August 7.742 UT), and lNelson et al. I d2007) observed a sim- 
ilar supersoft X-ray phase on day 789 (2007 August 28 UT). 
To search for the end epoch of supersoft X-ray phase, i.e., 
fo-hiiming, we use the automatic analyzer in the Swift web page 
(Evan s et al. 1120091) . The X-ray (0.3-10 keV) count rate with 
the Swift XRT are plotted in Figure[24] We have checked that 
the early several data points of the automatic anal yzer are very 
similar to those obtained bv lNess et al. I d2007d) . The epoch 
of ?H-buming is somewhere between 820 and 1070 days after 
the outburst, being ro ughly consistent with the prediction by 
lHachisu & Katol d2007ft . 

iSala et alTl (12008) found, from their X-ray spectrum fit, that 
oxygen/neon-rich WD atmosphere models provide a better fit 
than carbon/oxygen-rich WD atmosphere models. Therefore, 
we adopt here the chemical composition of "Ne nova 3" and 
recalculated model light curves. Our new best-reproducing 
model is a 1.07 Mq WD (Figure [24b. from which we obtain 
Wak = 131, f w ind = 355, f H -burnin g = 1000 days as listed in Ta- 
ble^ The estimated WD mass of 1.07 Mq is very close or 
nearly equal to the upper limit ma ss for CO WDs bo rn in a 
binary, i.e., M co < 107 Mq (e.g. JUmeda et al. Ill999t) . being 
consistent with the suggestion that the WD of V5116 Sgr is 
an ONeMg white dwarf. 

Adopting a magnitude of /«i/, max = 7.9 at maximum (to = 
3.6 days) along our model light curve of the 1.07 Mq WD, 
we obtain "intrinsic" ti = 15.6 and t^ = 28.5 days as listed in 
Table [6] The V magnitudes are roughly fitted with our model 
light curve at least until day ~ 100. 

We have checked the accuracy of WD mass determination. 
We obtain a best-reproducing model of Mwd = 1.0 Mq for 
X = 0.55 ("Ne nova 2"), and M WD = 0.85 M Q for X = 0.35 
("CO nova 2"). Then, we have 



Mwd(X) « {1.07 + 0.7(X-0.65)} M, 



0! 



(74) 



for V5 1 16 Cyg. Only when the hydrogen content is X ~ 0.65, 
the WD has a mass as large as 1.07 Mq, which should be 
consistent with an ONeMg WD. On the other hand, when the 
hydrogen content is as low as X ~ 0.35-0.55, then the WD is 
less massive (0.85-1. 0M Q ) and probably a CO WD. 

From the FF method, we obtain the distance modulus to 
V5116Sgr,i.e., 

[(m-M) v ] FF = m w -M w = 18.2-2.0= 16.2, (75) 

for the 1 .0 Mq WD ("Ne nova 2") model. Two other attendant 
models of 1.05 Mq and 0.95 Mq WDs ("Ne nova 2") give 
similar results of (m-M) v = 17.6-1.6= 16.0 for 1.05 M© or 
(m-M)v = 18.7-2.3 = 16.4 for 0.95 M Q WD. 
The LC method gives 



[51og(£/)+A v 



I V1974 Cyg 



[5\og(d)+A, 



I V5116 Sgr 



-3.9. 



(76) 

With (m-V)v = 12.3 from Equation (O for V1974 Cyg, we 
obtain 

[{m-M) v ] hC = [5log{d)+A v ] V5U6 Sgr = 16.2. (77) 

Kaler- Schmidt's law with £3 = 28.5 days and Delia Valle & 
Livio's law with t>i = 15.6 days, and both with my. max = 7.9 at 
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FIG . 25 . — Same as Figure [22l but for V574 Pup. We plot three model light 
curves of Mwd = 1.1 {thick dash-dotted lines), 1.05 (thick solid lines), and 
1 .0 Mq (thin solid lines) WDs for the envelope chemical composition of "Ne 
nova 3." We select the 1.05 Mq WD as a best model among three of 1.1, 1.05, 
and 1 .0 Mq WDs and can reproduce the supersoft X-ray data. Optical data 
are taken from AAVSO, Variable Star Observers League in Japan (VSOLJ), 
and IAU Circulars 8443 and 8445. The X-r ay flux data (large diamonds) 
with error bars (magenta bars) are taken from Ness et al. 1 2007a) and from 
an automatic analyzer of the Swift web page I Evans et al. 20091). 



maximum (fo = 3.6 days) yield 



[(m-M) 



V.maxJ MMRD 



= [51og(d)+A v 



I MMRD 



= 16.0(16.2), 
(78) 

being consistent with the FF and LC results. Then the distance 
modulus is (m-M)v = 16.2 ±0.2 and the distance t o V5116 
Sgr is d = 12± 1 kpc for A v = 0.81 dBurlak 1 120081) . These 
di stance estimates ar e consistent with d = 1 1 ± 3 kpc discussed 
bv lSalaetal.l(l2008l) . 

The o rbital period of Pn r b = .1238 days (2.97 hr) was de- 
rived by Dobrot ka et al. I (120081) from the orbital modulations 
with a full amplitude of 0.3-0.4 mag. Then, we estimate the 
epoch when the companion emerges from the nova envelope. 
From Equation d49b , we obtain M2 = 0.25 Mq, which yields 
a = 1.15 R Q and R^ = 0.58 R for M x = 1.07 Mq and, as a 
result, Emerge ~ 91 days at /? p h ~ a. 

6.9. V574 Pup 2004 

V574 Pup (Nova Puppis 2004) was discovered indepen- 
dently by Tago and Sakurai on 2004 N ovember 20.672 UT 
(JD 2453330.172) at mag about 7.6 dNakano et al. 1 12004b . 
This object was not det ected on November 12 and 16 with 
limiting mag about 13.6 (Nak ano et al. 112004b . so we assume 
that t 0B =JD 2453326.0 (November 16.5 UT) is the outburst 
day. The optical and near infrared l i ght cu rves are plotted in 
Figure [25] Although ISiviero et al. I ( 120051) derived m Vm . dx « 
8.0 on JD 2453335.6 (November 26.1 UT), this is not'the V 
maximum but the second peak as can be seen in Figure [25] 
We derive m v , max « 7.0 at t w JD 2453332.0 (November 22.5 
UT), 6 days after th e outburst based on the data reported in 
IAU Circular 8445 dSostero et al. II2004I) . After the V mag- 
nitude drops by about 2 mag at t « JD 2453333.7 (Novem- 
ber 24.2 UT), it increases to 8.0 mag at t w JD 2453335.6 
(November 26.1 UT). After that, the V magnitude declines to 
about 9.0 mag and stays for 20 days and then declines almost 
monotonically. In this way, V574 Pup shows an oscillatory 
feature in the very early phase. 

Swift XRT observations of V574 Pup showed that the 
nov a became a supers oft X-ray source before 2005 August 
17 dNess et al. Il2007al) . that is, 274 days after the outburst. 
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iNess etal.1 d2007bl) reported that V574 Pup was still in a su- 
persoft X-ray phase on 2007 December 7 and 13, that is, 1121 
days after the outburst. To search for the end epoch of super- 
soft X-ray phase, we use the automatic analyzer in the Swift 
web page (lEvans etal. 1120091) . The X-ray (0.3-10 keV) count 
rate with the Swift XRT are plotted in Figure [25] We added 
the last observational point from Swift web page. The epoch 
of fH-buming is somewhere between 1130 and 1320 days after 
th e outburst. 

iRudy et al~l (|2006) suggested that the observational proper- 
ties of V574 Pup were intermediate between those typical of 
explosions on a CO or an ONeMg WD. Their Spitzer obser- 
vation shows that V574 Pup reveals the strong coronal lines 
of Mg and Ne ions (but no Ne II at 12.8 microns), the behav- 
ior of which is suggestive of an explosion on the surface of 
an ONeMg WD. Therefore, we have calculated light curves 
for three different chemical compositions, i.e., "Ne nova 3," 
"Ne nova 2," and "CO nova 2." However, only the case of 
"Ne nova 3" provides a reasonable fit to the observation. Our 
best-reproducing light curve is plotted in Figure [25] for the 
chemical composition of "Ne nova 3." The WD mass is esti- 
mated to be ~ 1.05 ±0.05 Mq, which suggests a possibility 
that the WD is an ONeMg WD, because the estimated WD 
mass of 1 .05 M Q is marginally consistent with the theoreti- 
cally obtained upper li mit mass for CO WD s born in binaries, 
Mco < 1.07 M (e.g. JUmeda eTaT1[T999l) . 

For the other two chemical compositions, i.e., "Ne nova 2," 
and "CO nova 2," we do not find any reasonable light curves 
which simultaneously satisfy both the start and end of super- 
soft X-ray phase, that is, the duration of supersoft X-ray phase 
is too short to be compatible with the observation, like as has 
already shown in Figure l20l 

We determine tz and tj, along our model light curve. As- 
suming that our model light curve reached its maximum at the 
discovery day (fo = 4.172 days after the outburst), we obtain 
"intrinsic" 1 2 = 18.6 and tj, = 35 days. The other timescales are 
obtained to be f bre ak = 144, f wind = 390, to-burning = 1 140 days 
(Table©. 

From the FF method, we obtain the distance modulus to 
V574 Pup, i.e., 

[(m-M) v ] FF = m w -M vl = 17.4-2.0= 15.4, (79) 

for the 1 .0 M© WD ("Ne nova 2") model, although the mod- 
els of "Ne nova 2" are not shown in the figure. Two other 
attendant models of 1.04 M Q and 0.95 M© WDs ("Ne nova 
2") give similar results of (m-M)v = 16.9- 1.7 = 15.2 for 
1.04 M Q or (m-M)v = 17.8-2.3 = 15.5 for 0.95 M Q WD. 
The LC method gives 



[51og(£//10)+A v 



I VI 974 Cyg 



[51og(d/10)+A v 



I V574 Pup 



= -3.2. 
(80) 

With (m-V)v = 12.3 from Equation dT5} for V1974 Cyg, we 
have 

[(m-M) v ] hC = [51og( fl '/10)+A v ] V574 pup = 15.5. (81) 

Kaler- Schmidt's law with t^ = 35.0 days or Delia Valle & 
Livio's law with tz = 18.6 days and my max = 7.6 at maximum 
(fo w 4.2 days) give 



[(m-M) V[ 



= [51og(d/10)+A v ] 



MMRD 



= 15.5(15.7), 
(82) 

being consistent with the FF and LC results. Thus the distance 
modulus is (m-M)y = 15.5 ±0.2 and the distance to V574 
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FIG. 26.— Same as Figure [22] but for V458 Vul. Optical V (filled red 
squares), visual (small red open circles), R (ocher circles with plus), and / 
(magenta open diamonds) are taken from AAVSO and VSOLJ. The X-ray 
flux data (filled star marks for hard and open large diamonds f or sup ersoft 
X-ray) with error bars (green bars) are taken from Ness et al. (2009). We 
plot three model light curves of Mwd = 0.95 (thick dash-dotted lines), 0.93 
(thick solid lines), and 0.9 Mq (thin solid lines) WDs for the chemical com- 
position of "CO nova 4." We obtain f wl nd — 370, 416, and 456 days, for the 
0.95, 0.93, and 0.9 Mq WDs, respectively. Therefore we select the 0.93 Mq 
WD as a best model among three of 0.95, 0.93, and 0.9 Mq WDs by fitting 
our model value of ? w ; n( i with the observed turn-on of supersoft X-ray, i.e., 
'wind ~ 'x-on = 415 ± 15 days, because supersoft X-ray is self-absorbed by 
wind itself. We also indicate the epoch when the companion emerges from 
the photosphere and the break point of optical light curve in order to explain 
the flux of hard X-ray (filled star marks). 
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FIG. 27. — Same as Figuref26]but for chemical composition of "Ne nova 
2." We plot three model light curves of Mwd = 10 (thick dash-dotted lines), 
0.95 (thick solid lines), and 0.9 Mq (thin solid lines) WDs. We obtain f WIn d = 
352, 430, and 531 days, for the 1.0, 0.95, and 0.9 Mq WDs, respectively. 
Therefore we select the 0.95 Mq WD as a best model among three of 1.0, 
0.95, and 0.9 M WDs. 

Pup is d = 4.6 ± 0.3 kpc for A v = 2.2 dBurlakll2008l) . 

6.10. V458 Vul 2007 

V458 Vul (Nova Vulpeculae 2007) was discovered by Abe 
on 2007 August 8 54 UT (JD 2454321.04) at mag about 
9.5 dNakano et al. Il2007al) . This object was not detected on 
July 23, 31, and Au gust 4 with limiting mag about 11.5 
dNakano et al. ll2007cA . so we assume that t 0B =JD 2454317.0 
(August 4.5 UT) is the outburst day. The light curve is shown 
in Figure [26] The V magnitudes show a striking multi-peak 
(or multi-pulse) structure in the early phase. This kind of 
multi-peak structures were observed in HR Del 1967 and 
V723 Cas 1995, both of which are considered to be a nova 
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on a CO WD. So we assume here that V458 Vul was also a 
n ova on a CO WD. 

iDrake et al. I ((2007) reported the supersoft X-ray detection 
of V458 Vul with the Swift X RT on 2007 October 18 .05 UT 
(75 days after the outburst). Tsu iimoto et al. I (12008), how- 
ever, concluded, with the Suzaku observations, that the spec- 
trum can be described by a single temperature (~ 0.64 keV) 
thin thermal pla sma model, suggest ing shock-origin of the X- 
ray. Recently, Drak e et aL~l y008) reported, with the Swift 
XRT observations obtained from 2008 June 18 to Septem- 
ber 16 (about 320-410 days after the outburst), that V458 
Vul have entered a new phase characterized by a highly 
variable supersoft X-ray component accompanied by par- 
tially anti-correlated variations in the ultraviolet. They also 
pointed that an earl ier report on entry into the supersof t phase 
(IDrake et ai~ll2007l) has proven premature. iNess et al. 1 (2009) 
reported a summary of the early phase X-ray observation with 
Swift in which a supersoft X-ray phase started about 400 days 
after the ou tburst. Figur e [261 shows the supersoft X-ray fluxes 
reported in Nes s et al71 (120091) . The spectra show hard {star 
marks) components similar to Tsujimoto et al.'s until day 
~ 400 and then turned into supersoft {open diamonds). 

The chemical composition of V458 Vul is not yet reported. 
So we here assu me "CO nova 4 " in Table[T| whi ch is similar t o 
that of HR Del dTvlenda 1119781) and V723 Cas dlijirnall2006l) . 
Our model light curves are plotted in Figure [26] This nova 
has an oscillation phase in the early V, R, and / light curves. 
We thus fit our universal decline law with the bottom line of 
each oscillation as shown in the figure. The bottom line of 
/ magnitudes is nicely fitted with our model light curve until 
day ~ 100. On the other hand, the other two, V and R magni- 
tudes, are gradually departing from the model light curve due 
to emission line contributions. 

In the previous subsections, we have determined the WD 
masses mainly from X-ray light curve fitting because optical 
and near IR model light curves (i.e., free-free emission model 
light curves) have a similar shape of light curves and are not 
sensitive to the WD mass as explained in Figure [8] In the 
case of V458 Vul, only the emergence of supersoft X-ray was 
clearly detected but the supersoft X-ray phase does not end 
yet. Therefore, we have to estimate the WD mass only from 
the fx-on time. 

The supersoft X-ray flux replaced the hard X-ray compo- 
nent and rose up around t ~ 415 days as shown in Figure [26] 
As shown in V382 Vel, we proposed a hypothesis that hard 
X-ray component is a shock-origin by optically thick winds. 
Therefore, the substantial flux of hard X-ray indicates that 
optically thick winds still blow during the hard X-ray phase. 
On the other hand, supersoft X-ray may be self-absorbed by 
wind itself. Therefore, the emergence of supersoft X-ray indi- 
cates that the op tically thick wind has st opped or substantially 
weakened (e.g., Hachisu & Kato 2003). 

Our supersoft X-ray model light curves show a gradual rise 
in Figure [26] because our flux is the photospheric flux and 
does not include the effect of absorption by winds. Here we 
suppose that supersoft X-ray is absorbed during the optically 
thick wind phase. The sudden rise in the supersoft X-ray flux 
just after the optically thick winds stopped matches our view 
of de nse optically thick winds (see, e.g., Hachisu & Kato 
120031 for a sharp increase/decrease supersoft X-ray flux of 
RXJ0513.9-6951). 

The date of supersoft X-ray emergence is sharply deter- 
mined to be f wind = 415 ± 15 days for V458 Vul. We have 
calculated three models of 0.95, 0.93, and 0.9 M© WDs and 



obtained f w i nd = 370, 416, and 456 days, respectively. Then 
the WD mass is estimated to be M WD ~ 0.93 ± 0.03 M Q 
for "CO nova 4" from the emergence time of supersoft X- 
ray, i.e., f w ind ~ ?x-on ~ 415 days. Here we regard that the 
sharp increase in the supersoft X-ray flux corresponds to the 
end epoch of optically thick winds at about day 415. For 
the best model of M WD = 0.93 Mq among 0.95, 0.93, and 
0.9 M© WDs, we obtain f brea k ~ 145 days, f wind ps 415 days, 
and fH-burnin" ~ 1250 days (2011 January 5.5 UT) as listed in 
Table [6] 

We have checked the accuracy of WD mass determination 
especially for CNO abundance. We adopt other chemical 
composition of "Ne nova 2," same hydrogen content X = 0.55 
but different CNO abundance Xcno = 0.1 vs. 0.2, and we 
obtain f wind = 352, 430, and 531 days for M WD = 1.0, 0.95, 
and 0.9 Mq WDs, respectively. Here we adopt the best- 
reproducing model of Mwd = 0.95 Mq among these three WD 
masses as shown in Figure [27] This result may suggest that 

Mwd(Xcno) « {0.93 -0.2(X CNO - 0.2)} M Q , (83) 

for V458 Vul. The slope of -0.2 for CNO abundance is 
much smaller than that (~ 0.6) for hydrogen content. There- 
fore, if we regard the usual range of chemical composition 
as Xcno = 0.1-0.3 and X = 0.45-0.65, the ambiguity in 
WD mass determination (±0.02 Mq) for CNO abundance 
is much smaller than that (±0.06 Mq) for hydrogen con- 
tent. Therefore it is likely that M WD = 0.93 ± 0.08 M Q for 
a broad range of the chemical composition Xcno = 0.1-0.3 
(X CNO = 0.2 ± 0. 1) and X = 0.45 - 0.65 {X = 0.55 ±0.1). 

We have estimated tj and tj, along our model light curve. 
Placing our model light curve on the bottom line of multi- 
pulse light curves as shown in Figure [26] and assuming that 
the V magnitude attains its maximum of ~ 9.4 mag at to ~ 4.4 
days, we obtain "intrinsic" ti = 20.7 and £3 = 38.7 days. The 
assumed value of our theoretical peak is ~ 1.4 mag dimmer 
than that of the observed peak (8.0 mag). 

From the FF method, we obtain the distance modulus to 
V458 Vul, i.e., 

[(m-M) v ] FF = m w -M w = 19.3-2.3= 17.0, (84) 

for the 0.95 M Q WD ("Ne nova 2") as shown in Figure [27] 
Two other attendant models of 1 .0 Mq and 0.9 Mq WDs give 
similar results of {m-M) v = 18.8-2.0= 16.8 for 1.0 Mq or 
{m—M)v = 19.9 — 2.7 = 17.2 for 0.9 Mq WD. Here we do not 
use the model light curves of "CO nova 4," because they are 
not calibrated yet. 
The LC method gives 

[5 logtf/10) + A v ] V19?4 Cyg - [5 logtf/10) + A v ] V458 Vul = -4.7 

(85) 

With (m-V)y = 12.3 from Equation for VI 974 Cyg, we 
have 

[{m-M) v ] LC = [51og(rf/10)+A v ] V458 m = 17.0. (86) 

Kaler-Schmidt's law with £3 = 37.2 days or Delia Valle & 
Livio's law with t-i = 20.4 days, and both with my timK = 9.3 at 
maximum (to = 4.4 days) yield 

[(^-MV, max ] MMRD = [51og(<//10)+A,] MMRD = 17.1 (17.3), 

(87) 

being consistent with the FF and LC results. Thus we obtain 
the distance modulus of {m-M)y = 17.1 ±0.2 and the distance 
of d = 11 ±1 kpc to V458 Vul for Ay = 3AE(B-V) = 3.1 x 
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0,6= 1.86 dLvnch et al. ll2008tlPoggianill2008l:IWesson et al. I 
2008). This distance is roughly co nsistent with d = 13 ± 3 kpc 
estimated bv lWesson et al. I ( 120081) . 

Finally, we introduce some characteristic features concern- 
ing binary nature, that is, t he emergence of the companion 
from the WD photosphere. iGoranskij et al. I (120081) reported 
the orbital period of P or b = 0.58946 day with a full ampli- 
tude of 0.4 mag. The orbital light curve shows a saw-tooth 
shape similar to those in other supersoft X-ray novae, e.g. 
CI Aql and V723 Cas. Because Equation d49b is not appli- 
cable to such a long orbital period of P or b = 14.15 hr, we 
adopt a condition that the companion mass is smaller than 
0.74 Mq. This condition comes from q = M 2 /M\ < 0.8 
(Mi = Mwd = 0.93M Q ) required for thermal ly stable mass 
transfer in a binary (e.g., Webbi n¥eFari[T983l) . If we assume 
M 2 = 0.74 Mq as an upper limit of the secondary mass, then 
we obtain a = 3.5 Rq and R* = 1 .4 Rq. It is f e merge ~ 62 days 
at R p h ~ a. Even if we take a smaller mass of M 2 = 0.5 Mq, 
we obtain similar values of a = 3.3 Rq, R* = 1.45 Rq, giving 



emerge 



■ 65 days at R^ ~ a. 



Tsuiimot o et al. I (|2008 ! ;|) explained the early Swift X-ray 
data as shock-origin, which started about 70 days after 
the outburst and peaked at about 140 days and then de- 
creased as can be seen in Figure [26] The sho ck may be 
formed by collision between t wo ejecta shells (e.g., Friediung 
119871: ICassatella et al. 1 120041) or by collision between nova 
winds and the comp anion star dHachisu & Kato II2005L 120061 : 
lHachisu et ai~ll2008l) . 

If these X-rays come from the shock between the optically 
thick WD wind and the companion, the emergence of hard X- 
rays should be coincident with the emerge of the companion 
from the WD photosphere. This is roughly consistent with 
Emerge ~ 65 days. On the other hand, the break of free-free 
emission light curve is caused by a more steeper decrease in 
the wind mass-loss rate after the break point. Therefore, the 
drop of X-ray flux at/near 140 days after the outburst is also 
consistent with our model of fb rea k ~ 145 days. The hard X- 
ray lasted at least until about f w ; n d s=s 410 days, that is, until 
when the optically thick wind stopped. 

It should be noted that, for the orbital period of V458 Vul, 
R. Wesson et al. (private communication) proposed another 
orbital period of 98.1 min. If it is the case, the orbital pe- 
riod of P OT b = .0681 day (1.635 hr) is much shorter than 
that obtained by Goransk ii et al. I (12008) and this period gives 
M 2 =0.12 Mq from Equation (09), a = 0.71 R Q , 



,R*=QA0Rq, 

and femerge ~ 136 days (at /? pn ~ a). Then, the emergence time 
of the companion star is too late, that is, t emerge ~ fbreak ~ 140 
days, and our model on the origin of hard X-ray may not be 
viable. 

6.11. Results on prediction formulae and f? time 

In the previous subsections, we have fitted our model light 
curves with those of ten classical novae, in which a supersoft 
X-ray phase was detected, and estimated various nova param- 
eters. Figure|28]shows a summary of our results, that is, fx-off> 
fx-on, and t 2 (from top to bottom) against tj, for the ten novae, 
excluding V59 8 Pup but instead includi ng VI 974 Cyg, which 
was studied by lHachisu & Kato I (120061) . 

This figure shows a very good agreement between the data 
of "intrinsic" (t 2 , h) and Equation d29l ). This is not surpris- 
ing because Equation d29b can be derived from the universal 
decl ine law with F\ oc f~ 175 as shown in Figures [3] and [6] (see 
also lHachisu & Kato 2006), and both t 2 and f? are calculated 
from our model light curve of the universal decline law. 
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FIG. 28. — Supersoft X-ray on/off and "intrinsic" h time are plotted against 
"intrinsic" f 3 time for 10 novae, V1974 Cyg 1992, V382 Vel 1999, V1494 
Aql 1999#2, V4743 Sgr 2002, V574 Pup 2004, V5 1 16 Sgr 2005, V1281 Sco 
2007, V2467 Cyg 2007, V458 Vul 2007, and V597 Pup 2007#1. From top to 
bottom, fx-off ; tx- o n , an d int rinsi c ti in units of day. Solid lines: central values 
of Equations 13 II . 1301 , and j29t . Upp er/lower dashed lines: upper/lower val- 
ues of Equations <3 1 1 . )30K and {29}, respectively from top to bottom. Each 
value of fx-off ~ fa-burning (open squares), f X - n ~ Avind (open triangles), in- 
trinsic t% (open circles), and fj is taken from Table|6] 

Figure|28]also shows reasonable agreement of the observed 
fx-on and fx-off against f3 with Equations ( f30b and Pit , re- 
spectively. Theoretically, fx- on (~ the epoch when optically 
thick winds stop) and fx-off (~ the epoch when hydrogen-shell 
burning stops) are uniquely determined if the WD mass and 
the chemical composition are given. However, t 2 and f3 times 
are theoretically not unique even if both the WD mass and 
the chemical composition are given. These timescales depend 
further on the initial envelope mass as already explained in 
Section[3] Therefore, we must obtain t 2 and f3 times from each 
outburst of various novae. Good agreement between our pre- 
diction formulae and the observation on fx-on and fx-off sug- 
gests that our nova models of optically thick winds and free- 
free model light curves are reasonable for the entire descrip- 
tion of nova evolutions and, at the same time, the scatter in the 
initial nova envelope mass is not so large from their average. 
In our modeling, we estimated their average relations, t-$ vs. 
fwind or f3 vs. fH-buming, from the model of V1668 Cyg. Such 
a relatively small scatter of the 10 novae from these relations 
strongly suggests that we reasonably predict the supersoft X- 
ray phases of novae from this diagram, i.e., from Equations 
d30l l and OTb if we can estimate accurate "intrinsic" t 2 or f3 
time from early optical light curves. 

7. CONCLUSIONS 

Our main conclusions consist of three parts: the first part is 
for basic properties of our model light curves which may be 
summarized as follows: 

1. We model optical and infrared light curves with free- 
free emission, and UV 1455 A and supersoft X-ray 
light curves with blackbody emission, based on the op- 
tically thick wind mode l, as done in our previous paper 
( Hach isu & Kato 1 120 06). The model parameters gov- 
erning the light curves are the WD mass (Mwd), chem- 
ical composition of the WD envelope (X, Y, Z, Xcno, 
and X Ne ), and the initial envelope mass (M env ,o). 
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Distance Moduli of Novae 
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a distance modulus estimated from direct fitting (FF) with the calibrated fr 


ee-free model light curve 









distance modulus estimated from two light curve (LC) fitting 
c distance modulus estimated from Schmidt-Kaler's law of Equation (38) or Delia Valle & Livio's law of Equation {39} in parenthesis, 
i.e., from the Maximum Magnitude -Rate of Decline (MMRD) 

reference for quoted values of Ay or E(B - V), where we assume that Ay = 3. 1E 03-V): lfl teadetal. 1 2008), 2 Shore et al. i 2003), 
3 )Vanlandingham et al. I {200%, 4{Russell et al. I t2007al). 5{Ness et al. 1 12008d), 6fl iiima & Esenoglu i 2003). 7 Mazuk et al~J 120071) . 
8 Burlak (2008), 9- Wesson et al. (2008). 10 Chochol et al. 1 1993), 1 1 Stickland et al. (1981), 12 Downes & Duerbeck i 2000 ), 13- 
Wuetal. il989) 

e error means a standard deviation of the four distance estimates, i.e., FF, LC, and two MMRDs 

f 



parenthesis means the distance estimated from other method, see text for details 



g the distance is taken from Downes & Duerbeck i 2000) 
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FIG. 29. — Photospheric radius against photospheric temperature forWDs 
with various masses only during the static phase (after the optically thick 
winds stop). Red solid lines denote those for 0.6, 0.8, and 1.0 Mq WDs with 
the chemical composition of "CO nova 2." Blue solid lines denote those for 
1.0 and 1.35 Mq WDs with the chemical composition of "Ne nova 2." Black 
thick solid line indicate an approximate formula of Equation 4B H . Its ±0.1 
upper/lower bounds are indicated by two black dashed lines. Short horizontal 
lines indicate the WD radius (i.e., the bottom of a nuclear burning region) of 
0.6, 0.8, 1.0 and 1.35 Mq WDs. See text in Appendix B for more details. 



The peak brightness of our model light curve depends 
also on the initial envelope mass (M env ,o) even if we 
fix the WD mass and the chemical composition of the 
envelope. This indicates that nova speed class, i.e., ?2 
and ?3 times, depends also on the initial envelope mass, 
in other words, the mass accretion rate to the WD. 

Our UV 1455 A model light curves for different WD 
masses are homologous, that is, they are almost over- 
lapped to each other if we stretch them in the direction 
of time (f) by a factor of l// s and properly normalize 
the flux. All the free-free emission model light curves 
for different WD masses are also homologous and over- 
lapped with each other if we rescale them as t 1 = t / f s in 
the time and F' v = f s F v in the flux (see Figures[3]and[6]), 
that is, they are overlapped to each other in the t'-m' v 
plane, where m' v = -2.51ogf^+ const. The more mas- 
sive the WD, the smaller the timescaling (or stretching) 
factor of / s . We call this "a universal decline law." 

Using the universal decline law, the apparent mag- 
nitudes of various WD mass models are restored by 
my = Wy + 2.51og/ s and their absolute magnitudes are 



calculated from My 



, + 2.5logf s -{m-M) v . There- 



fore, we can determine the absolute magnitudes of free- 
free emission model light curves by calibrating their 
absolute magnitudes with the known distance moduli 
(m-M)y of V1668 Cyg for CO novae and V1974 Cyg 
for neon novae. 



2. Our free-free emission model light curves reasonably 
well reproduce optical and near IR light curves of no- 
vae. The decay timescale of our free-free emission 
model light curves depends strongly on the WD mass 
but weakly on the chemical composition. So the nova 
speed class is mainly determined by the WD mass. 



6. Since the absolute magnitudes of our free-free emission 
model light curves are now calibrated (Tables|2]and[3]l, 
we can estimate the distance moduli of individual no- 
vae by directly fitting them with observation. We show 
that the distance moduli thus obtained for VI 500 Cyg 
are consistent with those calculated from the nebular 
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expansion parallax method. This is a new method for 
obtaining the distance modulus to a nova. 

7. It has been empirically proposed that the absolute V 
magnitude 15 days after optical maximum, My (15), is 
almost the same among various novae. We show that 
there is only a small scatter of 0. 1-0.3 mag for the abso- 
lute brightnesses of our model light curves 15 days after 
maximum (see Figures 0] and |7J. This results theoreti- 
cally support the above empirical relation on My(15). 

8. We found that, when the timescales of two nova light 
curves are similar, the absolute brightnesses of the two 
novae are almost the same during the period in which 
the universal decline law is applied. Therefore, we are 
also able to estimate the distance modulus of a nova by 
comparing it with the other nova the distance of which 
is known and the timescale of which is similar. We also 
show that the distance moduli thus obtained for GK Per 
are consistent with those calculated from the nebular 
expansion parallax method. This is another new method 
for obtaining the distance modulus to a nova. 

9. Our free-free model light curve reasonably follows the 
observed light curve when free-free emission domi- 
nates the continuum flux in novae. Even when the nova 
has a transition oscillation or a characteristic secondary 
peak, our free-free emission model light curve can be 
fitted with the bottom line of observed nova light curves 
(see Figures [10][TT] andHSJl. 

The Second part contains timescales of novae, our predic- 
tion formulae of supersoft X-ray on and off times, and our 
theoretical MMRD relations, which are derived from the uni- 
versal decline law: 

10. We derive several relations between the characteris- 
tic timescales of classical novae, i.e., fbreak, fwind, and 
?H-burning, which depend mainly on the WD mass and 
weakly on the chemical composition but almost inde- 
pendent of the initial envelope mass. 

1 1 . On the other hand, ?2 and tj, depend sharply on the ini- 
tial envelope mass. Adopting a typical (average) case 
of the initial envelope mass, we obtain "average" rela- 
tions between ?3 (or tz) time and the three timescales 
of fbreak, ?wind, and fH-buming- These proposed relations 
are reasonably fitted with published data of novae. We 
propose the relation between the "intrinsic" tj, time and 
fwind, i-e., Equation d30b . as a prediction formula for the 
start of a supersoft X-ray phase, and the relation be- 
tween the "intrinsic" tj, and fH-buming, i.e., Equation d3Tb . 
as a prediction formula for the end of it. Here we pro- 
pose the "intrinsic" tj, time which is calculated along 
our free-free emission model light curve fitted with the 
observation. 

12. We obtain a relation between the maximum magni- 
tude, ffly jlnal = my max + 2.51og/ s , and the rate of de- 
cline, log?3 = log(?3/ s ), through the stretching factor / s 
of the free-free emission model light curves, i.e., Equa- 
tion (1321 1. A more massive WD corresponds to a smaller 
stretching factor of f s while the initial envelope mass 
determines a set of (f3,'«(/ max ). Calibrating the distance 



moduli (m-M)v of VI 668 Cyg and VI 974 Cyg, we de- 
rive usual MMRD relations between the absolute max- 
imum magnitude My. max and log tj, by changing / s for a 
given set of (?3 , »?y. max ), i.e., Equation (l34l . 

13. Our MMRD relations theoretically confirm that the 
main parameter is the WD mass, Mwd- A more massive 
WD with a smaller / s corresponds to a faster nova, i.e., 
a shorter time. However, there is a second parameter, 
i.e., the initial envelope mass M env ,o, corresponding to 
{t^m'y max ). The larger the envelope mass, the brighter 
the maximum magnitude, and the shorter the rate of de- 
cline . The main parameter Mwd produces a main trend 
of MMRD relation and the second parameter M env o can 
reasonably explain the scatter around the main MMRD 
relation. 

14. Our MMRD relations are also in goo d agreement with 
Kaler-Schmidt's law dSchmidt 1 119571) and Delia Valle 
& Livio's (1995) law at their central lines (values). 

The last part concerns light-curve fitting of ten classical no- 
vae in which a supersoft X-ray phase was detected: 

15. The universal decline law reasonably reproduces the 
observed light curves of ten classical novae, in which 
a supersoft X-ray phase was detected. The models re- 
producing simultaneously the optical and supersoft X- 
ray observations are WDs with 1.28±0.04M Q (V598 
Pup), 1.23 ±0.05 M Q (V382 Vel), 1.15 ±0.06 M Q 
(V4743 Sgr), 1.13 ±0.06 M Q (V1281 Sco), 1.2 ± 
0.05 M (V597 Pup), 1.06 ±0.07 M Q (V1494 Aql), 
1.04±0.07M Q (V2467 Cyg), 1.07 ±0.07 M Q (V5116 
Sgr), 1.05 ±0.05 M (V574 Pup), 0.93 ±0.08 M Q 
(V458 Vul). The first nine WDs are probably ONeMg 
WDs (ONe novae) while the last one is probably a CO 
WD (CO nova). Our light curve fittings also suggest 
various physical values of novae, such as the envelope 
mass of M env o at optical maximum, the ejecta mass of 
Mvind, etc. as tabulated in Table [7] 

16. The newly proposed relationships of Equations (f30b 
and (l3~TT l are consistent with the emergence/decay epoch 
of the supersoft X-ray phase of these novae including 
V1974Cyg. 

17. We estimate distances to the ten classical novae by us- 
ing our two new methods and show that they are consis- 
tent with those derived from the MMRD relation based 
on the "intrinsic" ?3 time as tabulated in Table [8] 

18. We propose a mechanism of early hard X-ray emis- 
sion originated from the shock between the optically 
thick winds and the companion star. Our model pre- 
dicts that the hard X-rays emerge when the companion 
appears from the white dwarf photosphere and decay 
when the optically thick winds stop. Our estimates of 
emerging-time f em erge of the companion from the nova 
photosphere are consistent with the observed hard X- 
ray appearance in V382 Vel, V1494 Aql, and V458 Vul; 
which may support our idea of hard X-ray emission. 
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APPENDIX 

SCALING LAW OF BRIGHTNESS AND ABSOLUTE FLUX OF FREE-FREE EMISSION LIGHT CURVES 

In this appendix, we present a formulation of scaling law for free-free emission model light curves and derive Equation ©. 



We determine the flux of our free-free light curves against that of 0.95 M Q WD (? 
proportionality constant in Equation (]]) as 



)} 



Mi 



{M WD } 



(I) 



V1668 Cyg). Now we explicitly express the 



(Al) 



In free-free emission, F v is independent of frequency is, then C^ MwD ^ is also independent of frequency because each term in the 
square bracket of the right-hand-side is also independent of frequency. The proportionality constant C^ MwD ^ probably depends 
not only on the WD mass but also weakly on the chemical composition and other physical conditions. However, we here assume 
that C^ MwD ^ depends only on Mwd when we fix the chemical composition. 

Then, we obtained F' v , = f % F v from Equation (HJl after time-stretching (t' = t/ f s and v' = f s v). Since the free-free flux is 
independent of the frequency v, i.e., F' v , = F' v , this simply means 



(A2) 



(A3) 



(A4) 



F jo.95M } (f) = /sF jM WD } (f//s) _ 

Now, we define the apparent V magnitudes as 

m v (t)=-2.5logFj M ^(t) + H v , 

and, for 0.95 M Q WD (» V1668 Cyg) with / s w 1, we have 

m'^t)=-2.5\ gF^ 9SM ^(t) + H Vl 

where we introduce a constant Hy in order to fit m! y of 0.95 M Q WD model with the V observation of V1668 Cyg, and Hy is the 
same among various WD masses. Using Equation (1A2| |. we have 

<(r)=m v (*// s )-2.51og/ s . (A5) 

By the way, we have already defined the apparent V magnitudes of our free-free emission model light curves as 

-| {M WD } 

A I- 

m'yit/fs) ; 



-2.5 log 



M wind 
V ph^ph 



+ Ky 



(A6) 



-I (<//.) 



in Equation ©. Here Ky is a constant defined to fit m' v of 0.95 M Q WD model with the V observation of V1668 Cyg, and 
common among all the WD mass models (from Mwd = 0.55 to 1.2 M Q by 0.05 M Q step) for a chemical composition of "CO 
nova 2." We saw that all the light curves are almost overlapped to each other in Figure|2] This simply means that 

{0.95 M Q } 

(A7) 





{Mwd} 




M wind 




_ V ph /? ph_ 


(f/f.) 


_ V ph^ph_ 



-I (!) 

Substituting Equation (lAll i into Equation (1A21 > and using Equation dA7b . we have 

£.{0.95 M s } ^^(HMwd} (Ag) 

From Equation (IA7t together with Equations (IA41 > and iA6\ , we derive 

m' v (t/f s )-Ky wm^(f) + 2.51ogC {0 ' 95Me} -// v . (A9) 
Therefore, we may regard that niy{t) = m' v (t/f s ), or in other words, 

<(*//.) = m v (f// s )-2.51og/ s , (A10) 

if we define Ky as 

Ky = H v -2.5 log c {a95Me} . (All) 

Thus, we derive Equation (0. Therefore, we plotted all the brightnesses of novae in the t' = t // s vs. m' v = my-2.51og/ s diagram 
as shown in Figure [2] 

It should be noted that the UV 1455 A and supersoft X-ra y mo del li ght curves follow F' v , = f s F u after stretching but do not 
satisfy F v i = F v , so that these fluxes do not obey Equations (lA2b and (lA51 l. because the blackbody emissivity depends on the 
frequency v and, as a result, the spectrum changes after time-stretching (v 1 = f^v). Therefore, we have to normalize the UV 
fluxes by different factors to fit the peaks with the observational peak of VI 668 Cyg, as already mentioned in Section [3] 
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AN EMPIRICAL FORMULA OF WHITE DWARF PHOTOSPHERIC RADIUS 

The photospheric radius of a typical nova during a supersoft X-ray phase is usually much larger than the Chandrasekhar radius. 
It is about ten times larger (~ 6 x 10 9 cm) than the Chandrasekhar radius at the beginning of a supersoft X-ray phase and quickly 
shrinks with the photospheric temperature being increasing. Therefore, we here present a simple formula on the radius of WD 
photosphere during a supersoft X-ray phase, which may be used for estimating the X-ray emission area. The photospheric radius 
of the WD envelope in the supersoft X-ray phase is plotted in Figure [29] against the photospheric temperature for various WD 
masses with the chemical composition of "CO nova 2" and "Ne nova 2." The left edge of each line corresponds to the epoch when 
the wind stops (at f w ; n d » fx-on = when the supersoft X-ray emerges). As the photosphere shrinks, the photospheric temperature 
rises. The rightmost point of each line denotes the epoch when hydrogen nuclear burning stops. We see that the photospheric 
radius R p h is almost inversely proportional to a square of the photospheric temperature 7^, i.e., R ? h oc l/7p 2 h , because the total 
luminosity is almost constant, L p h = AnR^oT^ ~ const., in the supersoft X-ray phase except the very later phase. 

The radius of a typical fast classical nova (1 .0 M Q WD with a composition of "CO nova 2") can be approximated by the black 
solid line as 

logfl ph (cm) = -2.07 (log7p h (K)-5.7) +9.27, (Bl) 

for 5.45 < logTph (K) < 5.8. The dashed lines denote the upper and lower bound values of log/? p h ±0.1 which covers the most 
of the WD mass representative for classical novae, i.e., from 0.8 M with "CO nova 2" to 1.35 M with "Ne nova 2." The slow 
novae are not covered in the region, i.e., 0.6-0.7 M Q with "CO novae 2" has smaller photospheric radius as shown in Figure |29l 
This figure also shows the WD radius, i.e., th e bottom of a hydrogen burning zone: the Chandrasekhar radius for Mwd < 1 .33 M Q 
and values taken from lNomoto et al. I (119841) forM WD > 1.33 M Q . 
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